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SUMMARY 

This paper, based partly on lectures on radio astron- 
omy given at the Columbia University Summer Institute in 
Space Physics in 1964, reports on the observations of radio 
galaxies and quasi- stellar radio sources (quasars). The 
basic physical processes responsible for continuous radio, 
emission are studied. /j. 
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* 


by 

Wolfgang Priester* 
Johan Rosenberg 
Goddard Space Flight Center 


INTRODUCTION 

The most challenging research problem in the field of radio astronomy at the present time is 
the attempt to understand the physical processes which cause certain extragalactic objects to emit 
radio radiation, the total power of which is up to six orders of magnitude larger than the radio 
emission of our own galactic system or of other "normal" galaxies. There are two kinds of those 
extremely strong extragalactic radio emitters known: the radio galaxies and the quasi- stellar 
radio sources. The latter name is nowadays abbreviated for convenience to "quasar" following a 
proposal by Chiu (Reference 1). 

The radio sources called radio galaxies often show a double or multiple structure where the 
radio emission is observed to originate in large areas far outside of the optically observed galaxy. 
There are, however, also objects where a single area of radio emission is nicely centered around 
the optically observed galaxy. There are two basic problems in connection with the radio galaxies: 
first, to explain the observed brightness distribution, in particular in those cases where there are 
several or at least two centers of radio emission separated by distances which are much larger 
than the diameter of the optically observed parent- galaxy, and secondly, to find the physical process 
which is capable of providing the extremely large amount of total energy required for the production 
of the observed radio power. 

The second problem is also the main problem for the quasi- stellar radio sources. Here the 
difficulties are even increased due to the comparatively small size of the areas of radio and light 
emission and due to the evidence that the total power emitted from these objects in the optical 
range is in the order of two orders of magnitude larger than the total light emission from the 
brightest giant galaxies. Further the radio emission of the quasars is equal to that of the brightest 
radio galaxies. 

Let us say a few words on the terminology of radio sources. When John Hey in 1946, (Ref- 
erence 2) discovered a point source in the constellation Cygnus, he called it a "radio star." In 


*This report is partly based on lectures on radio astronomy given by W. Priester at the Columbia University Summer Institute in Space 
Physics, 1964. 

'National Academy of Sciences— National Research Council Senior Research Associate with the Goddard Institute for Space Studies; on 
Leave from Bonn University. 
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the early 1950 ! s about 100 radio sources were observed, but none of them identified with 
either a peculiar star or another optically known object. Then, in 1952, Walter Baade and 
Rudolph Minkowski identified the source in Cygnus and a few others with somewhat peculiar 
galaxies (Reference 3). Since the sources were distributed isotropically on the sky— if we 
disregard the few supernovae remnants in our own galaxy which are observable radio sources— 
it became apparent that the great majority or the radio sources were indeed galaxies with 
abnormally high radio emission and hence the term TT radio star" was replaced by "radio 
galaxy. " So far, more than 70 sources have been identified with galaxies, the radio emission 
of which is several orders of magnitude larger than the one from a normal galaxy like our 
own or M 31 (Andromeda nebula). The term "radio galaxy TT is now generally used for these 
strong radio emitters. 

The term "radio star," however, reappeared again at the 107th meeting of the American 
Astronomical Society in New York in 1960 when Allan Sandage announced that the radio source 
3C 48 (No. 48 of the 3rd Cambridge catalog of radio sources) was identified with a 16th magnitude 
"star" in Triangulum by the perfect coincidence between optical and radio position. The identifi- 
cation was backed by the extremely small diameter of the radio source, less than 1". On a 90 
minute exposure with the 200 inch telescope, the "star" is accompanied by a faintly luminous 
nebulosity measuring about 5" x 12", which, however, does by no means resemble a galaxy. The 
spectrum was considered to be very peculiar since no unequivocal identification with lines of 
known atoms or ions could be made immediately. By early 1963, five radio sources (3C 48, 147, 
196, 273, and 286) were identified with similar star-like objects. All these radio sources have 
extremely small angular diameters of the order of one second of arc. The term "radio star" was, 
however, given up again, and replaced by quasi-stellar radio source when Maarten Schmidt (Ref- 
erence 4) discovered that the observed optical emission line spectrum of the source 3C 273 is 
explainable if one accounts for a redshift of AxA = 0.158 corresponding to an apparent recession 
speed of 47,400 km/sec. This pointed strongly towards an extragalactic object at a great distance. 
There was therefore no good reason anymore to believe that the actual size of these objects was 
comparable to the size of a normal star. 

In this paper we shall report essentially on the observational results on radio galaxies and 
quasars. Since there are no generally accepted views for the physical interpretation on the origin, 
evolution and energy supply of both kinds of objects, and since the proposed theories are likely to 
be subject to changes or improvement, we decided to leave out a detailed account of the proposed 
theoretical interpretations for the energy source of the objects. Instead we include a chapter 
treating the basic physical processes which are responsible for the generation of continuous radio 
emission. Both the radio galaxies and the quasars show typical non-thermal radio spectra which 
are generally believed to originate from magneto-bremsstrahlung of relativistic electrons spi- 
raling in a magnetic field ("synchrotron emission"). For didactic purposes and comparison we 
treat also the origin of thermal spectra from Coulomb-bremsstrahlung in the chapter preceding 
the chapter on the synchrotron emission. 
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BASIC DEFINITIONS 


Intensity and Flux 

The radiation intensity l v (6 ,<p) in a radiation field is defined as the amount of radiant energy, 
passing into the direction specified by (d 9 <p) through a unit area A, perpendicular to this direction, 
during one time unit into a solid angle of one steradian within a frequency band of 1 c/s. There- 
fore in c.g.s. units the dimensions of intensity are: 


[l v (i 9, <£)] - erg cm 2 ster” 1 sec 1 (c/s) 1 . 

Although not explicitly mentioned, I v (O, <p) also will depend in general upon time and the position 
of A. In Figure 1, then, the energy dE v , passing through the area da into the solid angle d fi during 
the time dt with frequencies between v and v + dv, is 


dE v - I v (£, <p) cos 6 dcr dt dfi dv . 

The radiation flux S v is defined as the amount of energy passing through an area of 1 cm 2 in all 
directions during one time unit with a frequency bandwidth of 1 c/ s. Hence the relation between 
flux and intensity is given by 


S, - I I v (9, 0) cos 9 dfi 
•'o 

with the dimensions 

[S v ] = erg sec” 1 cm -2 (c/s)” 1 

The symbol O denotes the whole sphere over 
which the integration has to be extended. If one 
expresses the element dfi in spherical coordi- 
nates (9, 0 ) (see Figure 1), we have dfi = sin d 6 d 0 , 
and henceforth the flux is 

/•2rr p7T 

= J J l v ( 6 , <p ) sin 6 cos 6 d& d<^ . ^ 2 ) 

The average intensity of the radiation coming 
out of the unit area into one hemisphere denoted 



Figure 1 — Test area da in a spherical coordinate 
system (6,<p). The dotted area is the projection of d^ 
onto a plane perpendicular to the direction ( 9,<p ), 
where dfi denotes an element of solid angle. 
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H is given by 


I 


<p) cos 6 dH 


1 


cos 6 dfi 


(4) 


If the test area is situated in a field of isotropic radiation , i.e., if i v (0, cp) is independent of (£, 0 ), 
then S v = 0. If there is no radiation coming in and the radiation going out is isotropic, we have 


F 


V 


S„ 

I v<0.*> = - • 


(5) 


The total flux over all frequencies is given by 


s 


f 


S M dv 


( 6 ) 


Its c.g.s. units are [S] = erg cm 2 sec 1 . 

Brightness Temperature 

We intend to find a relation between the emitted intensity at a certain frequency and the 
physical quantities determining the nature of the source. For this purpose it is convenient to 
compare the observed intensity with the intensity emitted by a black body at the same frequency. 

A black body is defined by the property that it absorbs all incident radiation at all frequencies. 
According to Planck's Law a black body having a temperature T radiates with an intensity given 
by 


B v (v, T) 


2hv 3 


1 


e h ^A T _ ! 


erg sec 


ster 


(c/s)“ 


where 


Planck T s constant: h 
Boltzmann 1 s constant: k 
velocity of light: c 


6.6256 ± 0.00 0 5 X 10 27 erg sec 
1.38053 * 0.00006 x 10~ 16 erg deg' 1 
2.997925 ± 0.000002 X 10 10 cm sec" 1 . 


(7) 
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The radiation of a black body does not depend on its material but only on the frequency and the 
temperature. The total radiant flux of a black body is given by Stefan-Boltzmann’s Law 


S 


f 

•'n 


S„ djy 


crT 4 , 


(8) 


where the Stef an- Boltzmann’s constant is a = 5.6697 x 10 " 5 erg cm " 2 sec " 1 deg" 4 . 


In radio astronomy the temperatures of the sources will normally be over 100°K. The fre- 
quencies in use will in general be less than 3.10 4 Mc/s (corresponding to \ > 1 cm). Hence the 
quantity h^AT is generally <0.015. So we can safely approximate e hv/kT by the first two terms 
of the Taylor expansion; 


e hv/kT _ ! 


hv 

kT 


(9) 


With this we obtain the Rayleigh- Jeans approximation of Planck’s Law for hvAT « 1: 


B 


R.J. 

v 


2hv 3 kT 
c 2 


2v 2 kT 


2k T 


erg cm" 2 sec" 1 


ster" 


(c/s)' 


( 10 ) 


This approximation is very important in radio astronomy and will be used throughout these 
chapters. Note that the quantum theoretical factor has disappeared. This formula was actually 
already derived by Lord Rayleigh, using classical methods. Inserting appropriate values for k 
and c we find 


B v r -J- = 3,07 x 10~ 25 i/ 2 T erg sec" 1 ster" 1 (c/s)" 1 cm" 2 


( 10 a) 


or 


B v R ' j * = 3.07 x 10 28 ^ 2 T Watt m 2 ster" 1 (c/s) 1 


( 10 b) 


in both cases taking M = Mc/s and [T] = °K. 


A similar approximation exists for extremely high frequencies, i.e., for h^/kT » 1. Then 
e hwkT _ i ^ e hvAT ^ 3 ^^ pianck’s Law becomes 



e ~hiVkT 


(11) 


which is called Wien’s approximation . 
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The brightness temperature T b O) of a source for a certain frequency is defined as the tem- 
perature a black body must have to emit an equal intensity at that frequency, so 

s B ( v - T b) ■ (12) 

Here \ v is the experimentally determined intensity. In radio astronomy, using the Rayleigh-Jeans 
approximation we find 


C x v 

T b( v) " 2k • ~2 (13) 

The brightness temperature corresponds only to the actual temperature if the source emits as a 
black body at that frequency, i.e., if the source is isothermal and if it absorbs all incident radi- 
ation. In all other cases the brightness temperature of a source will be a function of the frequency. 

RADIATIVE TRANSFER 

We will discuss here only a simplified case of radiative transfer applicable to most radio 
astronomical problems. A general treatment of the radiative transfer problem is given by S. 
Chandrasekhar (Reference 5). We will calculate the intensity of the radiation emitted by a layer 
of matter of finite thickness. The surface of the layer is assumed to be perpendicular to the line 
of sight. 

The absorption coefficient k v is defined as the fractional decrease of the intensity, when the 
radiation passes through a layer of 1 cm.: 

1 V - *v ds • (14) 

k v is generally dependent upon the position of ds, and upon the frequency. It is [* v ] = cm' 1 . The 
emission coefficient is defined as the amount of energy emitted by a unit volume per unit solid 
angle per second and per frequency band width of 1 c/s. Its dimensions are accordingly [e v ] = erg 
sec*" 1 cm" 3 (c/s)~ 3 ster” 1 . 

With the notation as explained in Figure 2 we can write the equation of radiative equilibrium 

l v ( s + ds ) dco dt dcr dv - I v ( s ) da) dt do* dv = - k v • I v ( s ) ds dco dt dr dv + € v ds dcr dt dw dv 


or 


dI v ' - lv ds + € v ds * (15) 

where i v , Kp and depend upon and in general on the position of ds. 
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In our simplified consideration we assume 
that T is a function of s only. We further make 
the presumption that all the radiation arriving 
at s comes from regions close enough to s to 
have their temperatures approximated by T(s), 
that is, that changes of the temperature over 
the mean free path length of the radiation are 
small. This condition is called local thermal 
equilibrium (LTE). It is a very useful restric- 
tion, because then Kirchkoff's Law holds, say- 
ing that the ratio of the emission coefficient to 
the absorption coefficient is independent of the 
material and given by Planck’s function: 



^(s)- 0 


Figure 2 — Simplified radiative transfer problem. The 
observer is assumed in the direction perpendicular to the 
surface of the layer, where s is the geometrical, and r v 
the optical thickness. 


* v = [T(s)j • (16) 

This equation holds independent of the way in which the LTE is established, whether by radiation, 
convection or conduction. Note that [T(s)] is independent of the direction of the radiation. The 
optical thickness t v is defined by 


&T V = - K v ds . (17) 

Clearly we want r v to be zero at s = s 0 , assuming there is no absorption between the layer and 
the observer (see Figure 2). And for s < s 0 , we want (s) > 0. Hence 

r° 

MO - J M s • (18) 

Note that r v is a dimensionless quantity. Division of Equation 15 by Equation 17 and using Equa- 
tion 16 yields 


d MO MO - - 

dr v (s) 1.(0 k v (s) Ms) ®«. [T(0] • (19) 

Multiplying Equation 19 by e ^ (s) : 


-t (s) dI v ( S ) - T v (0 _ _t v (s > d r-r/ 

e &r y (s)-e M s ) ' “ e L T ( s )J 

and integrating between the limits t v - 0 and r v = r v (0) (i.e., s = s 0 and s = 0) we find 


( 20 ) 


r\(°) 

K (0)e TvW - l v (s„) = - B„ r v {r v )e Tv &r v 

vo 


( 21 ) 
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Assuming an isothermal layer we have B v [T(s)] = B v (T) and hence B v (T) independent of 
r v ( s ) . Inserting this in Equation 21 yields 


I, (0) e' V0) - I, (s 0 ) = B„ (T) (e' Tv<0> - l) 


( 22 ) 


Assuming that no radiation enters the layer from the rear side we then have the intensity of the 
radiation leaving the surface of an isothermal layer perpendicular to the surface: 


K ( s o) 




(23) 


If T „ (0) 

layer [t v » 


» 1 so that 
l) is 


e~ T * <0) 


« 


1, then the intensity for the emission from an optically thick 


I 


V 



(24) 


Hence the brightness temperature of the outgoing radiation is equal to the temperature of the layer 


T b = T . (25) 

If r v (0) is small [r v « l) we approximate e Tv(0) by 1 - r v (0) and obtain the intensity of an 
emission from an optically thin layer 


l v = T , (°) B v - r v(°)7r kT - (26) 

Accordingly the brightness temperature of the radiation is for t v ( 0 ) « 1 

T b (v) - T v (0) T . ( 27) 

From Equation 23 it also follows directly that in general 

T b O) = T(l-e' T * (0) ) . (28) 

In an isothermal, homogeneous layer of ionized gas k v is obviously independent of position and 
direction; thus we have 


r s o 

T v(°) = ds = s o 


(29) 
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and 


l * ( s o) 


2v 2 kT 


(l-e" SS °) 


(30) 


From Equation 30 it is clear that the spectrum of a radio source will be determined by the depend- 
ence of the absorption coefficient k v on the frequency v. 

There are two physical processes known which are generally believed to cause the radio 
emission with continuous spectra from celestial objects. One is a thermal process (Coulomb- 
bremsstrahlung), the other a non-thermal one (synchrotron emission or magneto-bremsstrahlung 
of relativistic electrons). 

(1) The Coulomb -bremsstrahlung, i.e., the radiation emitted by electrons during Coulomb 
collisions, experienced when moving through the plasma with thermal (non-relativistic) velocities. 

(2) The Magneto-bremsstrahlung , i.e., the radiation emitted by electrons which are acceler- 
ated in a magnetic field, while moving with non-relativistic (cyclotron emission) or relativistic 
velocities (synchrotron emission). 

Thus we will derive two possible forms of spectra that radio sources may have. We can then 
compare observational data with these theoretical curves and decide whether the source is thermal 
(Coulomb-bremsstrahlung) or non-thermal (Magneto-bremsstrahlung), where we are still left 
with the possibility that the observed spectrum corresponds to neither of the two theoretical ones. 
All extragalactic radio sources (galaxies, radio galaxies and quasars) so far show spectra which 
can be interpreted as synchrotron emission spectra. We shall treat this process in a later 
chapter. First, for didactic reasons, we shall derive the thermal radio spectra from ionized 
gas-clouds. 


THE ABSORPTION COEFFICIENT FOR COULOMB-BREMSSTRAHLUNG 
The Energy Emitted in a Free-Free Transition 

In a completely ionized interstellar gas cloud without a magnetic field the radio frequency 
radiation with a continuous spectrum is merely due to free-free transitions of electrons in the 
Coulomb field of the ions. This radiation is called thermal radiation because the electrons in- 
volved have thermal velocities (i.e., energies of about 1 ev). As the radio frequencies will be 
emitted largely by electrons undergoing small accelerations, we will restrict ourselves to electrons 
having small deflections while passing a nucleus. This means the hyperbolae of the electron orbits 
are to be flat, which puts a lower limit on the collision parameter p (Figure 3). As the path on a 
flat hyperbola is nearly a straight line, we will disregard the y-component of the acceleration. 


From classical electrodynamics (Reference 6), we know that the power passing through a 
unit area normal to the direction of the radiation is given by the Poynting vector at the position of 
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Figure 3 — Free-free transition of an electron at a nucleus 
(Ze). The path of an electron passing a positively 
charged nucleus with a collision parameter p large 
enough to approximate the hyperbolic path by a straight 
line. 


the area element. So 

If = 131 = £ ffixH] . ( 31 ) 

In the case of the accelerated electron, moving 
with non- relativistic speed, Equation 31 
becomes 

dQ 1 e 2 sin 2 6 x 2 0 

da - c 3 r 2 er g cm sec (32) 

if the y-component of the acceleration can be 
neglected. In Equation 32 the electron charge 
is e - 4.80 X 10“ 10 e.s.u.; 6 denotes the angle 
between the direction of the emission and the 
acceleration, r is the distance between the 
electron and the area element under consider- 
ation. In Equation 32 we inserted x for v, 
according to our approximation. 

The total energy emitted per unit time by 
the electron is 


Q = 




r 2 s i n 0 d@ d<f> 


(33) 


The area element da - r 2 sin&dQdcp in spherical coordinates. The axis {6 = 0) of this system is 
given by the momentary direction of v. In our approximation this is always the direction of the 
negative x-axis. The integration yields 


Q 


3 c s 


l 2 erg sec -1 


(34) 


The Frequency Distribution 

The frequency distribution of the emitted radiation can be obtained from a Fourier analysis. 
As illustrated in Figure 3 we define the time t with t = 0 at the moment of the closest approach 
between the electron and ion (x = - p, y = 0). Then x(t) is symmetric in t, i.e., x(+t) = x(-t) , 
so we only need to use the cosine terms in the Fourier analysis of x. Here we follow closely the 
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analogous derivation by Kramers (Reference 7) (1923) for the x-ray spectrum. A detailed account 
can be found in the books of Unsold (Reference 8) or Jackson (Reference 6). Formally the Fourier 
analysis can be expressed by 


x(t) 


f 

•'o 


C (cc) cos (at ) dt 


where 


(35) 


C(c) 


1 


77 



x(t) cos (cj t) dt 


(36) 


Herein is ^ = 2w, v denotes the frequency. As Q contains the second power of x, we use Parseval’s 
theorem to obtain the frequency spectrum of Q . We thus find 

2e 2 

Q* da; = Qv dv p 77 C 2 <>) do, t (37) 


From general arguments in the Fourier analysis of a pulse of radiation, it is obvious that there 
will be hardly any contribution from frequencies with periods smaller than the duration of the 
pulse. Taking the characteristic transition time of the electron t 1 - p/v [thus we assume x to be 
negligible before the electron has reached A (x = 0, y = - p) or passed B ( x = 0, y = + p) in 
Figure 3], we find for the contributing frequencies: 

1 _ y 

Qi< r 1 - p • (38) 


As we have negligible emission for |y| > p, we have - p/v < t < p/u ; so during the transition kt | < 1 
and we can approximate cos ux by 1. Since the acceleration according to Coulomb’s Law is 
v = Ze 2 /mi 2 , we have 

.. _ Ze 2 

x ‘ ^£2 cos^(t) (see Figure 3) . ( 39 ) 

Inserting \/l - cosi//(t)/p in x we obtain 


Ze 2 , 

— - cos 3 ^(t) . 
mp 


(40) 


From Equation 36 we find by using the fact that the contribution of x(t) is only appreciable for 
1 1 1 <p/v, i.e., where cos (cut) * 1; 


C(«) 


1 

7 7 



x( t ) dt 



j cos 3 t ) dt 

mp 


(41) 
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We have tg^(t) = v • t/p , and thus dt = p/v cos 2 f d^. Here we assume that v is approximately 
constant during |t| < t x . Then Equation 41 becomes 




CO) 


Ze 2 

77mp 2 


f+rr/2 

COS 3 \p 

J-rr/2 


P 

v cos 2 ^ 


dp 


Ze 2 ♦ 2 ^ 
wipv 


( 42 ) 


Using Equation 37 and remembering that Q(o>)% 0 for o;> l/t, and that C(^) corresponds to Equation 42 
for oo < l/tj we find 


Q„ d„ 


2e^ 

3c 3 


/Ze 2 2\ 2 

r yT7mpv ) 


do) 


16 e 6 Z 2 
3c 3 m 2 p 2 v 2 


dv for 


277p 


(43) 


and 


Q v dv - 0 for * > ^ * ( 44 ) 

Note that in Equation 43, Q v is independent of v . The Coulomb-bremsstrahlung spectrum of a 
single electron will therefore approximately look like a step function. 


The Radio Spectrum of an Ionized Gas Cloud 

To find the spectrum of the radiation emitted by the electrons in the interstellar gas, 
we have to sum over the appropriate collision parameters and over the Maxwellian distri- 
bution of the velocities, since these are the only variables appearing in the spectrum of one 
electron. 

Let N x be the number density of the ions, and N e the number density of the electrons. Then the 
number (nj) of Coulomb collisions, with collision parameters between p and p + dp, during one 
second, that an electron with velocity v will experience, is equal to the number of ions in a cylin- 
drical ring of length v and of respective inner and outer radii, p andp + dp. So rij = (v * 2n p • dp)x N. 
collisions sec " 1 per electron. The number of collisions per cm 3 involving electrons with veloci- 
ties between v and v + dv and collision parameters between p andp + dp, dN (v, p) is then 


dN(v, p) - (v * 277 p * dp) x N. x (N e f(v) dv) collisions cm 3 sec 1 , (45) 

where f(v) denotes the Maxwell velocity distribution function. The energy emitted by a unit 
volume during one second in the frequency band between v and v + dv is then, using Equations 43 
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and 45 and the definition of the emission-coefficient 


■ r f 

Jp, Jo 


Q v (p, v) dN( v, p) 


f 2 f 16e 6 Z 2 

J J o v 5 Ni N ' f(v) 277 pvdpdv 


Jo “ J pt * 


32rr e 6 Z 2 N. N_ f” 
3c 


(46) 


From the kinetic theory of gases we know that the mean value of l/v is 1/v = VWwkT ; thus we have 


8e6z3N , N e ,/ST 1 Pj , 3 

-In — erg sec 1 cm (c/s) 1 . 


3J 


3c 3 m 


tt yjpf p, 


(47) 


We now have to insert appropriate values for p 2 and p t . First we discuss the lower limit p t . 

We restricted ourselves in the entire discussion to flat hyperbolae, i.e., a small angle of 
deflection. We will take the upper limit of the deflection angle 8 equal to 77/2 radians. Then as 

tg 5/2 ' Ze 2 /p ■ l/mv 2 we find 


Pi 



Ze 2 
3kT * 


(48) 


Another restriction on p x is the de Broglie wavelength of the electron. It has to be smaller than 
2 77 P jj otherwise our classical derivation does not hold. Hence 


* - h 

Pi ” 277 m u 


(49) 


Comparing Equations 48 and 49 for some typical values inH n regions (ionized interstellar hydro- 
gen clouds, which are the prototypes of thermal sources), we find, taking T = 10 4 °K, v = 700 km sec" 1 , 
z = 1 


pj ~ 6 x 10 8 cm, pj* - 1.6 x 10 8 cm . 

Generally we will find p* < p l9 and we will therefore use p 1# 

Now we discuss the upper limit p 2 . In the discussion of the Fourier analysis of the spectrum 
of one electron we found that there will only be appreciable radiation with frequency v for 
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2nv < l/t l = v/p . This imposes an upper limit onp so 

V 

P2 ” 27 TV (50) 

Furthermore, as we only want to consider two body collisions we want p < r/2 (where r is the 
average distance between the ions); r is given by 47Tr 3 1^/3 = 1 , so 


1 /477 \ ' 

2 (t n J 


3. 2 N. 1/3 


Another restriction on p is given by the following. Around every ion we can picture a region in 
which its Coulomb field is undisturbed by the surrounding electrons and ions; outside this region 
however the other charges will "screen" its Coulomb field. An electron, passing this ion with a 
collision parameter p larger than the radius D of this "screening," clearly is not affected by the 
ion (Reference 9a and 9b). Taking the Debye length for the screening radius D = (kT/^e 2 N e ) 1/2 , 
we obtain a third upper limit 


P>” = 6 9 (l hT [ “' (52 

In the regions we are considering we have generally p 2 ** » p 2 * . Example: in a H II region: 

T = 10 4 °K and N.^ N 3 < 10 3 cm -3 . Hence p* ^ 0.3 cm, p**^20 cm. Furthermore we find that for 
the H n regions, p 2 < p 2 * . In the corona of the sun, however, the density is much higher, so that 
there p 2 * < p 2 ; in that case we have to use p 2 * . We thus find for the emission coefficient, using 
Equations 47, 48, and 50, 


(2kT) 3/2 3 
2rr Ze 2 ym v 


erg cm** 3 sec" 


The derivation of Equation 53 holds only for two body collisions. Up to now no satisfactory treat- 
ment has been given of multiple scattering. We use Kirchhoff T s Law (Equation 16) to derive the 
absorption coefficient k v : 

4 e 6 Z 2 ^ 1 . ( ( 2kT) 3/2 \ 

K » ' 3cm 2 • N i N a v 2 T 3/2 ln \4.22nZe 2 1 ^ V J (54) 

The numerical values given in Equation 54 are taken from Elwert (Reference 10) and Oster (Ref- 
erence 11). The logarithmic term in Equation 54 is not sensitive for variations in its constant. 
For example, if T = 10 4 °K, and v = 2000 Mc/s, a change in the constant by a factor 10 changes 
k v only by 25%. 
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We are now in a position to determine the radio spectrum of an ionized cloud (HU region) by 
inserting x v into Equations 24 and 25 , and neglecting the slow variation with v of the logarithmic 
term in k v . In the case that the optical thickness r v = k u s 0 « 1 we have the intensity X v independ- 
ent of the frequency and proportional to T~ 1/2 : 


l v ~ v° * T " I/2 for r v « 1 . 
For large optical thickness r v » 1 we have 


(55) 


I 


V 


v 2 • T 


(56) 


In this latter case the intensity corresponds to the intensity of a black body at the same temperature. 
As an example we calculate the surface intensity of the radio radiation emitted by an average H 
n region. The average corresponds to a list of 9 H II regions investigated by Sharpless and 
Osterbrook (Reference 12) and Oort (Reference 13): T = 10 4o K, N. = N e = 10 cm' 3 , s 0 =34 pc. 

With these figures the optical thickness becomes r v = s 0 = 1.4 x 10 15 v' 2 where the frequency 
i/ is in (c/s). The optical thickness is equal to 1 for a frequency of v = 37 Me/ s. For smaller 
frequencies we have emission from an optically thick layer, therefore the intensity is proportional 
to v 2 \ for larger frequencies than 37 Mc/s we have emission from an optically thin layer, where 
the intensity is nearly independent of the frequency. The shape of the spectrum is presented in 
Figure 4. 


THE SYNCHROTRON EMISSION 

The electromagnetic radiation from ac- 
celerated electrons with relativistic veloci- 
ties was treated theoretically by Schott in 
1912 (Reference 14). In the laboratory the 
radiation from electrons circling in a syn- 
chrotron was observed first in 1948. The 
importance of this process for the interpre- 
tation of the non-thermal radio radiation, 
observed from the galaxy and from the ma- 
jority of the discrete radio sources, was 
realized by Alfven and Herlofson (Ref- 
erence 15) and by Kiepenheuer (Reference 16) 
in 1950. It is now generally believed that the 
non-thermal radio radiation is produced by 
Magneto-bremsstrahlung (synchrotron emis- 
sion) from cosmic-ray electrons which spiral 



FREQUENCY (me) 


Figure 4— Thermal spectrum. The theoretical radio 
spectrum of a thermal source, with T - 10 4o K, 
Nj = N e = 10 cm' 3 , s 0 = 34 pc. 


15 



» 


in the magnetic fields in our galaxy. Ginzburg investigated in 1951 in an important paper (Ref- 
erence 17) the relation between the observed spectral intensity distribution of the galactic radio 
emission and the energy spectrum of cosmic-ray electrons. At that time electrons were not yet 
observed in the cosmic rays. Ginzburg found that if the number of electrons is of the order of one 
percent of the number of cosmic-ray protons of the same energy then the observed radio spectrum 
could easily be explained by the synchrotron emission process. His prediction was nicely con- 
firmed by the observations of Meyer and Vogt (Reference 18 and 19) and of Earl (Reference 20 
and 21) who found a ratio of 3 ± 1 percent between the electrons and protons in the cosmic rays 
for energies above 500 Mev. 

In the first section of this chapter we shall describe the radiation pattern of an accelerated 
electron in general. In the remaining part we then will derive the relation between the syn- 
chrotron spectrum of cosmic ray electrons and their energy distribution. 

The Angular Distribution of the Radiation Emitted by an Accelerated Electron 

In the non- relativistic case, i.e., v/c « 1 , the power emitted by an electron accelerated in a 
magnetic field has the same momentary angular distribution as in the case of Coulomb- 
bremsstrahlung (refer to the previous section): 


dtr 


Am 


V 0 . 

— sin 6 erg sec ster 


(57) 


where 6 is the angle between the direction of the acceleration and the considered direction of the 
radiation, dfi is the element of the solid angle. 

In the relativistic case v/c£ 1, the angular distribution depends on v andv, and also on the 
general relativistic effect of the transformation from the rest frame of reference of the electron 
to the observers frame. We will first discuss the dependence on v and v and will then treat the 
relativistic effect on the observer. 

From the didactic point of view it is preferable not to restrict the discussion on the accelera- 
tion of an electron in a magnetic field, but to treat two general cases of acceleration: (1) where 
the acceleration and the velocity have the same or opposite direction (linear accelerator): v || v , 
and (2) where the acceleration is perpendicular to the speed vector (synchrotron): v _|_ v. 

We shall emphasize here the points which are important in our context, following essentially 
the treatment of the subject as given by Panofsky and Phillips (Reference 22) and Jackson (Ref- 
erence 6). We shall use the following abbreviations throughout: 


^ P = 


E 

m 0 c 2 
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hence 


1 

1 - f. 3 2 


where E is the total energy, m 0 the rest mass of the electron, e the speed of light. 


We treat first the case: v 11 v : 

The power radiated per unit solid angle 
has the following pattern (Reference 6): 


dP e 2 v 2 sin 2 t> 

jo - 7 7 ere sec ster 

^ 4t7 c 3 (1-3 cos#) 5 


(58) 


Note that now 6 is defined as the angle between 
the velocity and the considered direction of the 
radiation. The direction of v and v which here 
are either parallel or anti-parallel is clearly 
an axis of symmetry. In the non-relativistic 
case we had a toroidal radiation pattern with 
a sin 2 6 cross section. Now in the relativistic 
case the toroidal pattern is distorted and bent 
over into the direction of the speed vector, 
while maintaining that the emission is zero in 
the axis of symmetry. Illustrations are given 
in Figures 5 and 6. From Equation 57 the 
angle of peak intensity is easily found: 


90 ° 60 ° 



Figure 5— Radiation of a linearly accelerated electron. 
The angular distribution of the electromagnetic radiative 
power emitted by a linearly accelerated electron with 
fl- v/c = .50. The logarithmic scale has arbitrary units. 


cos 6 


= W (Vi+ is^-i) 


(59) 


For /3-lwe have 


1 

e m »* " 2y ' (59a) 

The power in the direction # max increases very 
rapidly with increasing 7 . It is proportional 
to 7 8 . The total power, obtained by integrating 
( 2 ) over all solid angles is 


P = 



(60) 


90 * 60 * 30 * 



90 * 60 * 30 * 


Figure 6— Radiation of a linearly accelerated electron 
(the same as Figure 5 with J3 = v/c = .99). The loga- 
rithmic scale has arbitrary units, different from those of 
Figure 5. 
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More important for the application in radio astronomy is the second case: v_L v . This is the 
case when an electron moves in a magnetic field where the field vector is perpendicular to the 
speed vector. We define a coordinate system in the following way: the z-axis in the direction of 
the velocity vector, the x-axis in the direction of acceleration. The corresponding polar coordinate 
system (c , 0 ) is determined by 0 measured from the positive z-axis, 0 in the xy-plane with 0=0 
in the positive x-axis. Then the power per steradian in the frame of the electron is given by 


dP 

an 


477 C 3 


1 

(1-/3 co s0) 3 


sin 2 6 cos 2 0 
y 2 (1-/3 cos O') 2 _ 


where again 3 - v/c , y - E/m 0 c 2 . 

For y » 1, /3~ 1, we can approximate (5) by 

dP _ 2 e 2 v 2 fi 1 47 2 # 2 cos 2 0 

dO = F ’ ”^3“ y (l +7 2^2)3 “ (l +y 2^2)2 


( 61 ) 


(62) 


In this case, we do not have an axis of symmetry. Here we find the maximum power at 0 = 0, the 
direction of the speed vector. We have two directions where no emissions takes place: in the 
xz -plane (<p = 0) if cos 0 - /3 , i.e., 0 % 1 /y. 

In Figures 7 and 8 the angular distribution of the radiation for /3 = 0.99 and y - 1.4 is given 
in the xz -plane (0 = 0) and the yz-plane (0 = 90°). It is interesting to compare this with the non- 
relativistic case and to note the way the toroidal pattern with v as axis is deformed by bending it 
over in the direction of v . This results in side lobes in the xz -plane. In the xz -plane the main 
beam is confined within an angle of 2 /y . It is therefore reasonable to define the characteristic 


90 ° 60 ° 40 ° 30 ° 20 ° 



90 ° 60 ° 40 ° 30 ° 20 ° 


Figure 7— Synchrotron radiation. The angular distribution 
of the electromagnetic radiative power emitted by a cir- 
cularly accelerated electron with/3 = .99 in theplane0 = O, 
i.e., the plane defined by v and v. The power is given 
by the radius vector in an arbitrary logarithmic scale. 


90 ° 60 ° 50 ° 40 ° 30 ° 20 ° 



Figure 8— Synchrotron radiation (the same as Figure 7 but 
now wi th 0 = 90°). The plane is perpendicular to the 
one defined by v and v. The scale is the same as in 
Figure 7. 
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width of the main beam by 


m o c 2 1 

e 0 - E ~ y ' 


( 63 ) 


For y » x, the power at 0 - 0 varies as y 6 . The total power over all solid angles is again obtained 
by integrating Equation 61. We then find 


P 


2 e 2 v 2 / E ) 

3 c 3 \ m o C 7 



erg sec 


(64) 


Note that the total radiation power of a linearly accelerated electron varies with y 6 , whereas the 
power emitted by an electron in a circular orbit varies with y 4 . However if we write Equations 60 
and 64 in terms of the change in momentum, i.e., the applied force, and use the relativistic ex- 
pression for force, we find P = (2e 2 /3m 0 2 c)(dp/dt> 2 for a linearly accelerated electron and 
P = ( 2e 3 y 1 / 3m 0 2 c ) (dp/dt) 2 for an electron in a circular orbit. Table 1 gives some values of Q 
with the corresponding values of y and the characteristic beam width. 


Table 1 


Properties of Electrons and Protons Moving in a Magnetic Field. 


Particle 

Types 

Velocity 

(km/sec) 

/S 

7 

^ 1/2 

T 

E (ev) 

R l (km) 

(forH = 10 -5 Gauss) 

Thermal 

10 2 

1 

3000 


- 

200°K 

5 x 10 s 

0.57 

Electrons 

10 3 

1 

300 

-1 

- 

2 x 10 4o K 

5 x 10 s 

5.7 


1.5 x 10 5 

*5 

1.14 

50° 

— 

6 x 10 s 

10 3 


2.7 x 10 s 

*9 

2.3 

25° 

- 

1.18 x 10 6 

3.4 x 10 3 

Relativistic 

Electrons 


,99 

7.1 

8° 

— 

3.6 x 10 6 

11.9 x 10 3 



~1 

2000 

5 x 10‘ 4 rad 

— 

10 9 

3.4 x 10 6 = ~10“ 7 pc 


~c 


2 x 10 12 

5 x 10" 13 rad 

- 

10 18 

3 x 10 15 = ~100 pc 

Thermal 

Protons 

23 

; 

7.5 x 10' s 

-1 

- 

2 x 10 4o K 

9.4 x 10 8 

250 


1.5 x 10 5 

.5 

1.14 

50° 

— 

1.09 x 10 9 

1.8 x 10 6 

Relativistic 

Protons 

~c 

.99 

7.2 

8° 

— 

6.6 x 10 9 

22 x 10 6 


~c 


10 9 

10" 9 rad 

— 

10 18 

100 pc 
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The Larmor Circle 

If an electron enters a homogeneous magnetic field, the orbit will be a. spiral in the no n- 
relativistic case, the radius of the circle determining the spiral given by 


m 0 V c 

R l - eH, cm • 


( 65 ) 


^ is the component of H normal to v. In the relativistic case m 0 would be replaced by m = m 0 y (see 
Equation 69). 

The angular frequency called Larmor frequency is given by 

eHi 

^h.o “ m Q c rad sec 1 (66) 


In the non- relativistic case the electron radiates like a dipole, with a frequency determined by the 
Larmor frequency; this emission is called cyclotron emission or gyromagnetic emission . Its 
frequency is (in Me): 


v h,o — 2.8 x Hi Me , (67) 

where H x is in gauss. Cyclotron radiation is circularly or elliptically polarized in the direction 
normal to H^ To obtain radiation with radio frequencies, i.e., 10 - 10 4 Me, we would thus need 
magnetic fields of 1 - 10 4 gauss, which are highly improbable field strengths in interstellar space. 
In the sun T s atmosphere, these field strengths are possible but a more detailed treatment has 
shown that only the harmonics of ^ H0 can be observed, not the fundamental frequency. 

In the relativistic case we have 


"H./3 



eHx 


1 

— rad sec 


- 1 


( 68 ) 


So "h ,/3 "h.o* We ^ US ^ in<d that the angular frequency decreases compared to the non-relativistic 

case, even though the velocity is increasing; hence the radius increases. As 2 ttR l = v/v, we have 
the Larmor radius 


R 


L 


V 

co 


m 0 cv 

~£h 7~ ^cm . 


(69) 
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Note that in Equation 69, for y » 1 , v % c; hence 

m o cv E ^ _E_ . 

^ eHi * m Q c 3 " eH l ’ (70) 

so in that case R is independent of the mass of the particle. 

Table 1 gives a few values for the Larmor radii of different particles for different energies. 
The field- strength used is 10“ 5 gauss. This value is obtained by assuming that there exists 
equilibrium between the energy density of the magnetic field and of the turbulent motions of the 
clouds in the galaxies. 

The Duration of the Observed Pulse of Radiation 

Equation 61 states that for y » 1, the radiation is emitted mainly into a narrow cone pointing 
into the direction of the instantaneous speed of the electron. From Equation 63 we know that the 
characteristic angular width of the cone is equal to 

^0 = 7 ' ( 71 ) 

In order to sweep across a distant observer the cone has to move from A to B in Figure 9. In 
the reference frame of the electron this will take time At = R £ 0 /c, where we assume /3^l. In 
other words the duration of the pulse observed by the observer will last in the reference frame of 
the electron (using Equations 70 and 71), 


At 



m 0 c 


(72) 


A signal emitted at time t at R(t) is received 
at time t* by the observer at T (Figure 9) 

t * 


and as I r - §( t ) | 

= , 

dt " c | r - R( t ) | ‘ dt (74) 

Here [r - R( t )] j\ r - R( t ) | is the unit vector into 
the direction of the observer, and dR(t)/dt 
= v(t). Hence dt'/dt = l - (v/c) cos 9 . As At is 


= t + 


]T-R(t)| 


(73) 


= |/[ r ~ R( t )] 2 , 

1 [? ~ R( t )] dR( t ) 



Figure 9 — Pulse time. The relation between the charac- 
teristic beam width and the characteristic time during 
which it is observed at a distance ?. 
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very small we can approximate At' by 


* 



At' 

- ( 

t 

1 - 

* 

cos At 


(75) 

For 0 Q 

« 1, 

cos 0 

0 

yi - - v/c 

= P, 

so the characteristic time 

for the 

observer 

is 




At' - 

/ 

v*\ 


/ m o c2 \ 2 

m 0 c 


V- 

rf- 

At 

= \ E ) • 

eH x 

(76) 


The Spectral Distribution 

According to general properties of Fourier integrals, the spectrum will have appreciable 
intensity up to a critical frequency given approximately by the inverse of the characteristic time . 
We will define the critical angular frequency as 

1 


3 eH i / i 
2 % c Uo 


rad sec ” 1 


(77) 


or 


CO 


c 


3 

2 ^ 


H, 0 


3 

2 


°>h.p y 


(78) 


Since the spectrum in fact does not decrease abruptly to zero at a certain frequency, we shall find 
slightly different definitions for the critical frequency in the literature. Here we used Schwinger’s 
definition, which is most commonly used in radioastronomy [see, for example, Oort (Reference 13) and 
Jager and Wallis (Reference 23)]. Ginzburg, however, used as critical angular frequency a value which 
is only 2/3 of the one given inEquation77. Jackson’s value is twice as large as o> c in Equation 77. 

By Equation 67 we have for the critical frequency v c = 4.2 Hj. ■ y 2 Me, where H is in gauss. In 
order to obtain radiation in the radio frequency range (10 - 10 4 Me) from relativistic electrons y 
has to be between 500 and 15,000 if we assume a magnetic field of 10‘ 5 gauss. This range for y 
corresponds to energies from 2.5 x 10 8 to 8 x 10 10 ev. It is interesting to note what energies 
would be required for relativistic protons if one wants to obtain synchrotron radiation in the same 
radio frequency range from them. From Equation 77 one can easily see that the proton energies 
must be larger by (m p /m 0 ) 3/2 , where m p is the rest mass of the proton. Therefore the energy range 
required for protons would be from 2 x 10 13 to 6 x 10 15 ev. 

The frequency distribution of the total energy of the pulse as observed by the observer is 
obtained by Fourier analysis of the emitted power. We omit here the details and refer the reader 
interested in them to the paper by Schwinger (Reference 24). 
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The frequency distribution is given by 


Hl F fe“) er S sec 1 (c/s) 1 
m Q c \ c/ 


eH. 


where v c = ^ = r 2 , and 


-fe) ■ £ 


K 5/3 ( x > dx ' 


(79) 


(80) 


k 5/3 is a special case of the Airy integral (or modified Bessel function or modified Hankel function). 
Oort (Reference 13) gives a graph of the frequency distribution. 

G. Wallis (Reference 25) found a very convenient approximation of the Airy integral. It is 
sufficiently accurate for all applications in radioastronomy: 

/ V \ 0 - 3 -v/v 

F w = 1.78 • ) • e ' . (81) 

F and F w have the same maximum values at the same value of :/: c ; furthermore f F w d? = 0.993 

/•co *0 

Fdv. Actually the spectrum of one electron is not continuous, but a discrete spectrum with values 

Jo 

of P v different from 0 only in the harmonics of the Larmor frequency. However as y > 500 and 
= 3/2 ol h 0 y 2 the harmonics of 0 are so closely spaced in the region of o, c that we can regard 
P, as a continuous spectrum (Reference 24). 

The Spectrum of the Radio Frequency Radiation 
Emitted by the Galactic Cosmic Ray Electrons 

In analogy to the observed energy distribution of cosmic ray protons we assume a differential 
energy distribution for the relativistic electrons as given by 

N e (E)dE = A • E~ g dE cm - 3 , [N e ] = cm -3 erg" 1 , (82) 

where N e is the number density of cosmic ray electrons with energies between E and E + dE, and g 
is called the spectral index of the energy spectrum, A is a constant determining the absolute numbers. 
If we now assume that along the line of sight, the magnetic fields in general are distributed at 
random, then we can regard the synchrotron radiation as being isotropic, if we take an average 
over a large range. P v (E) is the total energy emitted per sec by an electron, so we can say that 
4/7 1^ (E) = P u (E) , where \ v (E) denotes the average intensity of the radiation of an electron in any 
direction, i.e., i v (E) is the isotropic intensity. We now want to calculate the intensity of the syn- 
chrotron radiation emitted along the line of sight as received by an observer. Note that the in- 
tensity of radiation is not dependent upon the distance between the observer and the source, as 
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long as we assume that there is no absorption along the line of sight. Jager and Wallis (Ref- 
erence 23) confirmed that absorption is negligible for the wavelengths and densities we are con- 
cerned with. Combining the spectrum emitted from a relativistic electron having the energy E 
with the energy distribution of the electrons N e , we obtain the intensity of the radiation by inte- 
grating over the line of sight and over the cosmic ray energy range: 

r r P * , 

I„ - ~ 4 ~ r N e (E, r ) dE dr erg cm 2 sec ster 1 (c/s) 1 . (83) 

Jo J E Q 


In order to estimate the synchrotron radiation from our galaxy for a certain direction we assume 
for simplicity that N e (E) is independent of r , for 0 < r < r 0 , where r 0 denotes the border of the 
galaxy along the line of sight, and N e (E) = 0 for r > r 0 . Then Equation 83 becomes 



FREQUENCY (me) 

Figure 10— Non-thermal spectra. A few examples 
of radio spectra of non-thermal radio sources (Ref- 
erence 26a). 


K = P„ ( E ) N e (E) dE . (84) 

\ 

If we observe a radio source which emits syn- 
chrotron radiation and appears under the solid 
angle n, then the received flux s v is 

s , = i„ n 

r r p * r 

“ n I 4 tt N e i E ( x ’ y. z )] dEdxdydz (85) 

Jvoi Je q 

Here we assumed already a random distribution 
of the magnetic fields in the source, and the ab- 
sence of absorption throughout the source, i.e., 
that the mean free path of the photons is much 
larger than the diameter of the source. Assuming 
again that N e (E) is independent of (x, y, z) we find 

00 

n voi I 

s , = ~ 4 rr P„(E)N e (E)dE , ( 86 ) 

Je o 

where vol designates the volume of the source. 

Returning to the intensity of the galactic syn- 
chrotron radiation Equation 84 we obtain by 
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using Equations 79 and 82 for the spectrum P v and the differential energy distribution N e , 
respectively: 


I 


V 



( 87 ) 


We now use the Wallis approximation for y[v/v q ) , and abbreviate /3 e 3 /m 0 c 2 with C. We 
furthermore abbreviate the expression for the critical frequency v c - eE 2 where e stands for 
(3/47T)^H/m 0 3 c s ). Then Equation 87 becomes 


1.78 

4 TT 


, CAHv 0 - 3 e' 0 * 3 E"^ +0 - 6 > e 


r 


v/ -€E 2 


dE 


( 88 ) 


If we substitute E = iv/e x we obtain 


1.78 

4 77 r o CAH 


(7) <E_ll/2 f x-<* +0 - 6 > e-V'x 2 dx 


(89) 


This reveals that the intensity of the synchrotron radiation should be proportional to v~ (g ~ l)//2 
provided that our assumptions are valid. In fact, the observed spectra of the non-thermal com- 
ponent of the galactic radio radiation and the non-thermal extragalactic sources can be reasonably 
well represented by a simple power law: 


(observed) - v a . (90) 

In order to prove that the non-thermal radio radiation stems from the synchrotron process, 
we have to show that if we equate the observed spectral index a to the theoretical value (g - l)/2, 
and also the observed intensity at one selected frequency to the theoretical intensity, we obtain 
reasonable values for both the spectral index of the energy spectrum (i.e., "g") and the number 
densities of the cosmic ray electrons. 

Figure 10 shows 7 typical spectra of non-thermal radio sources (Reference 26b). The spectral 
indices are (from top to bottom): a = 0.63, 0.69, 0.59, 0.71, 0.70, 0.90 and 0.88. The average 
spectral index taken over a larger sample of typical sources is 0.7. The spectrum of the galactic 
radio radiation is slightly less steep. From the survey by Baldwin (Reference 27) at 81 Mc/s, 
Droge and Priester (Reference 28) at 200 Mc/s and McGee, Stanley and Slee (Reference 29) at 
400 Mc/s, one finds a = 0.61. The same value was obtained by Westerhout (Reference 30a) from a 
comparison of a larger number of sky surveys.* If we take a = 0.7 we have to require that the 
index g of the energy distribution of cosmic ray electrons in the non-thermal radio sources is 


*In 1962 Turtle et al. (Reference 30b) again investigated the spectra of the radiation from the galactic halo and the galactic disc. They 
obtained a slight change of the spectral index a. with the wavelength: a smaller value at the longer wavelengths. For simplicity we 
shall only discuss a constant value for the spectral index. 
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g = 2.4. This value is in reasonably good agreement with the observed energy spectrum of 
cosmic ray protons in the energy range from 10 9 to 10 10 ev. In the literature the integrated 
energy spectrum N* of the cosmic ray particles for energies larger than a certain limit E 0 is 
usually given. It is 


(E 


■<.) ■ f». 

•/E 


(E) dE 


- r-G 

r. 


(91) 


where G = g - 1. 

Cosmic ray electrons were discovered in 1960 by Meyer and Vogt (References 18 and 19) with 
an electronic detector and by Earl with a cloud chamber in a sky-hook balloon at an altitude of 
about 36 km (Reference 21). The lower energy limit was E 0 = 5 X 10 8 ev. The observed flux was 
f = 32 ± 10 el/m 2 ster sec. If we assume isotropy, we are able to calculate N e * from 

N c * = ^ f - 1.3 x 1(T 12 cm' 3 (92) 


for E > E 0 = 5 X 10 8 ev. 

We now assume that the same power law, as observed for the energy distribution of cosmic 
ray protons, also applies for the electrons. Hence for 10 9 ev < E< 10 10 ev, we take g = 2.4. Using 
Equation 91 we find: 


N* (E> 5x 10 8 ev) 


— — F -1.4 

1.4 L o 


1 . 3 x 10“ 12 cm -3 


(93) 


or 


A = 8.2 x 10 17 erg 1 * 4 cm 3 . 

We are now in a position to calculate the intensity of the galactic synchrotron radiation for a 
selected direction and frequency and to compare the result with the observed intensity. For the 
value g = 2.4 the integral in Equation 89 is immediately calculated if we integrate from 0 to <». 
One can easily see that the error introduced by integrating from 0 can be neglected. 


c 00 c 00 

x -3 c _1 / x dx = 2 e -u du - • 

Jo Jo 

Inserting the constants into Equation 89 and remembering that e involved H, we find 

l v ~ 2.0 x 10 -26 H 1 ' 7 r Q v "°- 7 . 


(94) 


(95) 
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If we now take H = 10 5 gauss and r 0 = 2 x 10 4 parsec = 6 x 10 22 cm, corresponding to a galactic 
longitude l 11 = 30° we have 


I v = 5.8 x 10 ~ i8 (j~j^) er e cm" 2 ster 1 sec 1 (c/s )”* . (gg) 

Observations in the 200 Mc/s survey by Droge and Priester (Reference 28) revealed for l 11 = 30°, 
a brightness temperature T b = 500°K. This corresponds to an intensity of I v (observed) = 6 x 10" 18 
erg cm" 2 ster” 1 (c/s)' 1 sec" 1 . This is to compare with the calculated theoretical value of 
l v = 5.8 x 10” 18 erg cm” 2 ster" 1 (c/s)" 1 sec" 1 . The extremely good agreement between the 
observed and theoretical value is, of course, accidental, since rather crude assumptions went 
into the theoretical derivation. Nevertheless, it proved that the synchrotron process is able to 
explain the observed spectral distribution of the galactic and extragalactic non-thermal radio 
radiation and also that the observed absolute intensities are reasonably well represented under 
certain assumptions for the energy distribution of the cosmic ray electrons and for a galactic 
magnetic field of 10” 5 gauss. 

The Lifetime of Relativistic Electrons 

The energy losses of fast electrons moving through interstellar matter (H atoms and ions, 

N = 0.1 cm" 3 ) are due to the following processes: 

a) magneto-bremsstrahlung (cyclotron and synchrotron emission), 

b) bremsstrahlung due to Coulomb collisions between electrons and protons, 

c) ionization of atoms, 

d) the inverse Compton effect, i.e., scattering of fast electrons on thermal photons of the 
insterstellar radiation field in the galaxy. 

For electrons with energies E < 10 10 ev, loss due to d) « loss due to b). For electrons with 
energies E > 10 10 ev, loss due to the inverse Compton effect is comparable to the loss due to 
Coulomb bremsstrahlung, but still small compared to the loss due to magneto-bremsstrahlung in 
a magnetic field of 10” 5 gauss, which is believed to be characteristic for the interstellar space. 
Because the atom densities we are concerned with are very small, and the electron energies very 
large we can neglect the loss due to ionization compared with the loss due to magneto- 
bremsstrahlung. Table 2 [from Ginzburg (Reference 17)] gives a comparison between the energy 
losses due to the respective processes. 

a) the lifetime of an electron , calculated on the basis of energy loss by synchrotron emission. 
The total power emitted is given by Equation 64: 


P 



erg sec 


(64) 
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Table 2 


Energy Losses of Electrons in ev sec 1 due to the Following Processes.* 


Electron Energy in ev 

Ionization 
N = 0.1 cm -3 

Coulomb 
Bremsstrahlung 
N = 0.1 cm" 3 

Synchrotron 
Radiation 
H = 10“ 5 Gauss 

10 8 

2.8 x 10~ 8 

8 x 10' 9 

4 x 10" 9 

10 9 

3.3 x 10" 8 

8 x 10" 8 

4 x 10" 7 

5 x 10 9 

3.7 x 10~ 8 

4 x 10“ 7 

10" 5 

10 10 

3.9 x 10‘ 8 

8 x 10' 7 

4 x 10“ 5 

5 x 10 10 

4.3 x 10‘ 8 

4 x 10" 6 

1 

O 

rH 


* After Ginzburg (Reference 17). 


Remembering that v = v 2 /R = v • u = v (eH/m 0 c) (m 0 c 2 / E), we find for v = c. 


P = 6.2 X 10~ 27 X H 2 X E 2 erg sec -1 , 

where E is in ev and H in gauss. So 


(97) 


dE 

dt ” 


- P 


3.8 x 10“ 6 H 2 E 2 Bev sec" 1 


(98) 


where E is in Bev andH again in gauss. Integration of Equation 98 between t = 0 and t = t , let- 
ting the initial energy be E 0 yields: 


1 

E(t) 


3.8 x 10 6 H 2 x t + r — 
E o 


(99) 


Hence 


E(t) 


3.8 xlO _6 H 2 E 0 t+ 1 


erg . 


( 100 ) 


The half-life of the electron being defined by E(t 1/2 ) a 1/2 E 0 , we find 


- 1/2 


3. 8 x 10" 6 x H 2 E. 


8. 4x 10 

H 2 x E„ 


-2 


years 


(101) 


where E is in Bev and H in gauss. 

Table 3 gives the electron half-lives for different values of H and the initial 
energy E 0 . 
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Table 3 


Electron Half-Lives and Critical Frequencies for Several Electron Energies. 




Electron Energy 
in ev 

1 1/2 in Years 

v in Me/ s 

c 

K 

c 

h = icr 5 

Gauss 

10 8 

8 x 10 9 

1.6 

^200 m 



10 9 

8 x 10 8 


^2 m 



10 10 

8 x 10 7 

1.6 x 10 4 

cm 



10 12 

8 x 10 5 

1.6 x 10 8 

~2 x 10 4 A 



10 14 

8 x 10 3 

1.6 x 10 u 

o 

-2 A 

H = 1(T 3 

Gauss | 

10 13 

8 

1.6 x 10 12 

-2 A 


RADIOGALAXIES AND QUASI-STELLAR RADIO SOURCES 
General Description 

Although the first detection of radio waves from outer space was in 1932 by Jansky, it was 
not before 1946 that a discrete source was observed (Reference 2). The first tentative identifi- 
cations with optically known objects were made by Bolton, Stanley and Slee (Reference 31) who 
identified the radio sources Taurus A, Centaurus A and Virgo A with the Crab nebula (NGC 1952, 
NGC 5128 and NGC 4486 respectively). But it was only in 1952, when Baade and Minkowski identi- 
fied the radio source Cyg A with a peculiar galaxy which seemed to consist of two galaxies in 
collision and also studied several other objects extensively, that astronomers generally accepted 
the identifications (Reference 3). Up to the present time there are 1159 sources listed by Mills, 
Slee and Hill (Reference 32) and 328 in the revised 3rd Cambridge catalogue (Reference 33), of 
which so far only 70 have been identified with visible objects. We will discuss in this chapter only 
the extragalactic radio sources that have been identified with optically observed objects. Many 
radio sources will probably never be identified because their distances are too large for the 
associated object to be found on the photographic plates. 

Normal galaxies like our own and the Andromeda Nebula (M 31) have a total power in the 
radio frequencies (between 10 and 10 4 Me) between 10 38 and 1G 40 erg sec” 1 . Their radio lumi- 
nosity is comparatively low. They are therefore difficult to observe when they are at large dis- 
tances. Of about 50 identified sources only 8 were found to be normal. The term ’’normal” is 
here defined to mean that the ratio of optical to radio emission is in the same range as for the 
closest neighbor galaxies. It is expected that the majority of all galaxies will be weak radio 
sources, i.e., normal galaxies, and that only a rather small group constitutes the ’’radio galaxies” 
which are much easier to detect due to their extensive radio luminosity. The total power emitted 
from radio galaxies in the frequency range from 10 to 10 4 Me extends up to 10 45 erg sec They 
therefore constitute the majority of the observed radio sources. 
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It is now generally accepted that the main process which generates radio emission with con- . 
tinuous spectra in the meter wavelength range in galaxies is synchrotron emission. Matthews, 
Morgan and Schmidt (Reference 34) investigated how identified radio sources are distributed 
among the different types of galaxies. They studied 52 objects. Their results are described in 
more detail in Table 4. 

The weak sources are mostly spiral and irregular galaxies like our galaxy, the Andromeda 
nebula, and the Magellanic Clouds. They are so far always single radio sources with a simple 
structure, centered in the galaxy, but often surrounded by an extensive halo. 

The strong sources so far seem to be all related to certain special types of galaxies, namely 
to the D galaxies, DE type galaxies, dumbbell galaxies, N galaxies and to the quasi-stellar objects 
which have much smaller dimensions than the bright galaxies. 

a) D type galaxies have the following characteristics: a bright nucleus, mostly not flattened, 
which is surrounded by an extensive envelope. The very large giant D galaxies are often the 
dominating objects of the clusters of galaxies. D galaxies frequently have diameters three to four 
times as large (40 - 80 kpc) as the other members. Their optical luminosities cover a range of 
a factor 10, but their radio luminosities differ within a range of 10 4 . Although all have the dis- 
tinctive nucleus and the extended envelope, there are large differences in size, the nuclei some- 
times consist of two or more condensations, sometimes the cluster of galaxies itself is elongated 
in the same way as the envelope of the dominating D galaxy. Often absorbing dust bands are 
visible, revealing rotational symmetry. However, no optical feature indicates whether or not a 
giant D type galaxy should be a strong radio source or not. 


Table 4 

Range of Radio Luminosity for Different Types of Galaxies.* 


Galaxy Type 

Number 

Percentage 

% 

Radio Luminosity 
in erg sec" 1 

Optical Power 

Weak Sources 





Spirals, Irregulars 

8 

15 

10 38 - 10 40 

average galaxy 

Ellipticals 

2 

4 

10 40 - 10 41 

10 42 erg sec" 1 

Strong Sources I 





DE type 

13 

25 

10 40 - 10 43 


D type 

16 

30 

10 42 - 10 45 

giant galaxy 

db type 

5 

10 

10 41 - 7 - 10 43 

10 44 erg sec" 1 

N type 

4 

8 

10 43 


Quasars 

4 

8 

10 44 - 10 45 

10 46 erg sec" -1 


*The classification is made according to Morgan’s scheme (Reference 35). 
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The correlation between D galaxies that are radio sources and those which belong to clusters 
is interesting. 50 percent of all D type radio sources belong to clusters, but only 5 percent of all 
D type galaxies in clusters are radio sources. This might perhaps indicate that a D type galaxy 
becomes a strong radio emitter only for a rather short period, possibly because of some catas- 
trophe occurring inside its nucleus. If we assume the lifetime of a D type galaxy (or better its 
nucleus as it characterizes the galaxy) to be 10 9 years as a reasonable estimate, then we find that 
it would spend 5 x 10 7 years as a radio source. 

The weakest of the D, DE type galaxies show a simple radio structure, the more luminous ones 
are identified with double or more complex sources. 

The diameters of the D sources (or their components) are, just like the case of the weak 
sources, comparable to the optical size of the parent galaxy. Only very large halo sources and 
very strong double or multiple sources extend much farther; the separations of the components 
of the double sources are up to about ten times larger than the size of the galaxy. 

b) The DE type galaxies indicate a transition between the elliptical and the D galaxies, i.e., 
these include all variations of the nucleus slowly fading into the envelope. It is remarkable that 
galaxies of DE type, closely resembling the D galaxies have a higher absolute radio luminosity 
than the ones resembling the elliptical galaxies. 

c) dh type galaxies (i.e., dumbbell galaxies) are very similar to the D type, the only difference 
is that there are two about equal very distinct nuclei contained in the same envelope. Their radio 
luminosity range is about equal to that of D sources. Three of the radio sources identified with 

db galaxies are single and only two are double radio sources. This indicates a rather large pro- 
portion of single sources compared with the general case of radio galaxies where only 5 to 10 
percent are single radio sources. However the statistical sample is still rather small. 

d) The N type galaxies have very bright and small nuclei, the radiation whereof constitutes 
most of the optical luminosity. They have very faint nebulous envelopes, often just visible as a 
streak, and never far extended. Their radio luminosity is on the average higher than the lumi- 
nosity of the D sources, however this may be due to selection effects. They are closely related 
to the compact galaxies observed by Zwicky (References 36, 37 and 38), characterized by a very 
small nucleus, only to be distinguished from stars by the lack of the diffraction pattern on the 
photographic plate. These compact galaxies seem to exist over a wide range of optical spectral 
types. Here too are no visible differences between the strong radio sources among them and the 
normal ones. Zwicky (Reference 38) suggested that the quasi-stellar radio sources are the most 
luminous, extreme cases of these compact galaxies. 

The N sources probably have a complex radio structure, with several small components 
which seem to be rather long and narrow (Reference 34). Since only a few have been found, 
and their angular diameters are very small, more refined instruments are necessary to 
resolve their structure well. 
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e) Quasi-stellar radio sources ("quasars") . The main optical features of the quasi-stellar 
radio sources are: 


1) they can only be distinguished from stars by their unusual excess brightness in the blue of 
the spectrum extending into the ultraviolet, 

2) their spectra show emission lines with similarities to the emission lines of planetary 
nebula. Often the lines are very broad and have a very large redshift, 

3) occasionally very faint wisps or jets can be optically distinguished on the photographic plates, 

4) their radio diameters are usually very small, about 1". 

We shall return to a more detailed descrip- 
tion of the quasars later in this chapter. 


Description of Individual Objects 

We will now discuss several of the more 
closely studied radio galaxies. 

Centaurus A, NGC 5128 

The galaxy NGC 5128, identified with the 
radio source Centaurus A, is a DE type galaxy 
with a very extended envelope. The distance of 
Cen A is about 4.7 Mpc. The nucleus is ellipti- 
cal with a broad absorbing band containing dust 
and gas in the equatorial region. This band is 
rotating about the major axis of the nucleus, like 
the nucleus itself, but with a larger velocity 
(Reference 39). 

The weak extensions of the optical features 
along the axis of rotation are very peculiar; they 
reach out to more than 40 kpc and turn into the 
region of the radio emission. The only plausible 
explanation for the extension of matter along the 
rotational axis seems to be the ejection of ma- 
terial (ionized gas) along a magnetic field that 
might be aligned with the rotational axis. 

The radio source shows a rather complex 
and very extended elongated structure. It ex- 
tends over about 9 degrees on the sky. This cor- 
responds to a diameter of 650 kpc. (The dia- 
meter of the Milky Way is about 30 kpc.) The 
detailed brightness distribution of Figure 11 
shows at least 3 double sources in the extended 



RIGHT ASCENSION (a, 962 ) 


Figure 11— Centaurus A (outer sources). The brightness 
distribution in Centaurus A at 21.3 cm. The small black 
dot in the center equals approximately the galaxy as 
shown in Figure 12 (Reference 42). 
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source, and the so-called M inner source” close to the galaxy. This might indicate that these com- 
ponents were ejected by the galaxy during 4 successive events. This galaxy represents a special 
class of radio galaxies— another outstanding member of this class is Fornax A, associated with the 
galaxy NGC 1316. These objects generally seem to have very large linear diameters in both their 
optical and their radio size. Their radio luminosities are in the range of 10 42 erg sec' 1 . Matthews, 
Morgan and Schmidt (Reference 34) give the following values 

Cen A: L = 10 4187 erg sec' 1 
For A: L = 10 41 * 27 erg sec" 1 

for the range 10 to 10 5 Me. 

The polarization of the radio emission has been carefully studied, together with the depolari- 
zation effects occurring in the outer envelopes of Cen A or along the line of sight (References 40, 

41 and 42). The inner source is highly polarized, in some wavelengths up to 13 percent, but the 
polarization is generally restricted to rather small areas. The two components are not polarized 
to the same degree as can be seen from Figure 12. Comparing the position angles of the axis of 
the nucleus, the double source closest to the nucleus and the magnetic field, we find differences, 
but they are confined within an angle of about 30°. 



Figure 12— Centaurus A (inner source). The central part of the radio source CentaurusA at a wavelength of 10 cm. 
The two drawings show the measured brightness distribution for two different polarization angles. In the left part the 
position angle of the electric vector is 25°, in the right part 115° (adapted from Reference 40). 
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Fornax A, NGC 1316 (D type galaxy) 

This source was identified by Mills (Reference 32). It is in many ways similar to Cen A. It 
also shows some absorption features across the nucleus that curve up on one side and down on the 
other. Arp (Reference 43) was able to obtain photographical evidence of very faint luminous ex- 
tensions. Further away from the galaxy they curve towards the centers of the double radio source 
connected with it. The diameter of NGC 1316 as it is usually seen on photographic plates (see 
insert, Figure 13) is of the order of 20 kpc, in the case of Cen A (NGC 5128) it is about 25 kpc. 

Polarization experiments showed that the magnetic field was not aligned with the line con- 
necting the centers of the two radio components. It is now clear that the magnetic field follows 
quite well the shape of the optical extensions. The linear size of the radio emission area of 
Fornax A can be expected to cover up to 200 X 300 kpc. Thus it is comparable to Cen A. The 
main difference in the radio structure is that Cen A has several double sources and the inner 
source which is very close to the nucleus, while For A consists apparently only of one extended 
double source (Reference 44) (See Figure 13). 

Perseus A, 3C84, NGC 1275 

The center of a strong radio source was found to be in the Perseus galaxy cluster and identi- 
fied with the peculiar galaxy NGC 1275. It appears a tightly wound spiral but it has an extended 
envelope and a very bright nucleus. Several condensations and outward radiating streaks are 

visible indicating violent events inside the 
nucleus of this galaxy. Its nucleus has typical 
"Seyfert" characteristics and shows strong 
emission lines exhibiting irregular features; 
the hydrogen Balmer lines being emitted from a 
region just outside the nucleus are very broad. 
This suggests that inside the nucleus large 
turbulent motions exist, and that hydrogen clouds 
are moving away from the center at high speeds 
of the order of 3000 km sec' 1 . As the galaxy is 
situated in a cluster, these high velocities might 
be due to a collision between two galaxies; 
however no nucleus of a possible other galaxy 
colliding with NGC 1275 can be seen. 

The radio source Per A has a very small 
core, somewhat smaller than the galaxy (20 
kpc) and a very extended halo of about 90 
kpc or perhaps even much larger (Ref- 
erences 45 and 46). The central source ac- 
counts for 80-85 percent of the radio lumi- 
nosities, the halo provides the remaining 
15 percent. 



Figure 13 — Fornax A (NGC 1316). The radio contours of 
Fornax A are from Wade (Reference 44). The large opti- 
ca! features are sketched from a photograph published 
by Arp (Reference 43). The square in the center contains 
the bright part of the galaxy NGC 1316. This part is 
given in the lower left corner ten times enlarged. There 
the brightness contours and the absorption patterns are 
sketched from the Cape Photographic Atlas of Southern 
Galaxies. The large optical features of NGC 1316 on 
the photograph envelope the SO p galaxy NGC 1317, 
the position of which is indicated by a dot 6.5' north of 
the center of NGC 1316. 
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A radio galaxy quite similar to NGC 1275 is NGC 1086 , associated with the source 3C71 . It is 
also a "Seyfert galaxy/’ its nucleus is very small and extremely bright, and the emission lines 
are broadened, indicating strong large scale movements in the center. The extended structure 
shows spiral arms, and there is some evidence of a hydrogen flare, escaping from the nucleus. 

The interesting feature of NGC 1086 is that it is a comparatively weak radio galaxy of ex- 
ceedingly small diameter, its radio luminosity is about 10 40 erg sec -1 , placing it between the 
normal and the radio galaxies, the diameter of the source (600 pc or less) is much less than the 
diameter of the galaxy (^15 kpc). 

Cygnus A is the brightest extragalactic radio source in the sky in the meter wavelength range. 
It was discovered by Hey (Reference 2). In 1952 Baade and Minkowski (Reference 3) identified 
this source with a faint object (apparent visual magnitude 15 ra ). The redshift in the emission lines 
is 0.056, corresponding to a recession speed of 16,800 km/sec. Using a Hubble constant of 100 km 
sec -2 Mpc" 1 , we obtain a distance of about 170 Mpc. The optical diameter is about 8 kpc. The 
optical spectrum of Cyg A is very similar to that of the Seyfert galaxies, although the broadening 
of the lines is not as large as in the case of NGC 1275 and the emission comes from a region 
larger than the emitting regions in other Seyfert galaxies. Its optical power in the emission lines 
is 2x 10 44 erg sec" 1 , its total power about 4x 10 44 erg sec' 1 (Reference 47). 

The radio source is double with two about equal components approximately 34 kpc in diameter 
and with a separation of 79 kpc. The components are extended along the major axis, on the center 
of which the galaxy is located. A polarization up to 8 percent has been measured. The direction 
of the magnetic field, however, is uncertain because of depolarization effects and Faraday rotation 
inside and outside the source (Reference 48). The radio luminosity is 10 44 7 erg sec" 1 . Cyg A 
is one of the brightest galaxies observed so far, in both the optical and the radio frequency range. 

Virgo A t NGC 4486 (M 87) 

Bolton, Stanley and Slee (Reference 31) identified Virgo A with the galaxy M87 (= NGC 4486), 
a member of the Virgo cluster. Baade and Minkowski (Reference 3) investigated M 87 on 
the Mt. Palomar 200 inch plates and especially its optical jet which protrudes from the 
nucleus. The galaxy is surrounded by a very large number of globular clusters. The nucleus 
of the galaxy is bright and small (-100 pc). The total envelope also encloses the jet and 
has a diameter of 6 kpc. The distance between the outermost condensation in the jet and 
the center is 1100 pc. As Virgo A is a strong source, polarization measurements were car- 
ried out on the jet to find out if this peculiarity is associated with the source. The jet condensa- 
tions were found to be optically highly polarized, but with slightly different directions (Ref- 
erences 49 and 50). The optical spectrum of the jet shows a blue continuum, without emission 
lines, similar to the Crab Nebula. This gives strong evidence that the optical radiation from the 
jet is also due to synchrotron radiation. There is also ionized gas present in the nucleus since 
there appear strong emission lines of OH that show double features, indicating that large amounts 
of gas are moving away from the nucleus at speeds that are, however, still comparable to the 
random motion inside the nucleus. 
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The radio source has a rather complex structure: it has been resolved in two components: 

1) an unpolarized core source which is double with a separation of 2.5 kpc (Reference 51). 

The core source is elongated in the direction of the jet, another indication of the importance of the 
jet as part of the radio source, 

2) a halo source emitting about 50 percent of the radio radiation. The halo is slightly polar- 
ized and has a diameter of about 20 kpc. 

3C33 

The peculiarity of this object is that it is a double radio source where the components have a 
very small size compared to their separation. It is a DE type galaxy with angular dimensions of 
3\ r 3 x 8V4. It has an outer envelope also of elliptical form (22" x 9"), but with its major axis 
aligned with the minor axis of the nucleus. The galaxy is probably rotating about the major axis 
of the envelope. The major axis of the double source forms an angle of 35° with the axis of the 
envelope. Polarization measurements of the radio source show that it is highly polarized, with 
the direction of the magnetic field in between the axes of the galaxy and the source. The separa- 
tion of the components is 200 kpc (3!8), the size of the components which are elongated along the 
axis is 7x14 kpc (8” x 16") (Reference 52). The strong polarization and the small size of the 
components show a highly aligned magnetic field which is strong enough to confine the relativistic 
electrons to a rather small region. 

Quasi-Stellar Radio Sources (Quasars) 

In the introduction we stated already that for several years (from 1960 through 1962) the 
quasi-stellar radio sources were believed to be peculiar stars within our own galaxy. This was 
due to their perfect star-like appearance on photographic plates. 

The main breakthrough in the identification came with Maarten Schmidt’s paper (Reference 4) 
on the radio source 3C 273. The source coincides in position with a rather bright (13 m ) star-like 
object. He was able to explain the optical spectrum when he accounted for a redshift of 0.158 
corresponding to a nominal velocity of 47,400 km/sec. After this discovery, 4 more sources be- 
longing to the same class were soon reinvestigated, as they had already been identified with 
stellar-like objects. These are 3C 48, 147, 196, 286. 

By December 1963 the number had increased to 9 objects. The additional 4 objects are 3C 47, 
245, 9 and 216. 3C 47, however, has a larger radio diameter equal to about 1 minute of arc.* The 
term "3C" refers to the 3rd Cambridge catalogue of radio sources (Reference 33). 

The identification of the source 3C 273 was facilitated by the fact that the spectrum shows the 
Balmer series down to H epsilon. The appearance of 5 Balmer lines is quite unusual for these 
kinds of objects. A further peculiarity is that the forbidden O II line at 3727 A, which is usually 
very strong, is extremely weak here. 

*At the second Texas Symposium on Relativistic Astrophysics (Dec. 1964) A. Sandage announced that 25 quasi-stellar objects have 
been found from their ultraviolet excess. 
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The star-like object is accompanied by a faint wisp or jet. The jet has a width of 1" - 2” and 
extends away from the star in position angle 43°. It is not visible within a distance of 11” from 
the "star” and ends abruptly at 20” from the "star.” 

Fortunately this radio source was occulted by the moon three times in 1962, on April 15, 
August 5, and October 26. The occupations were observed with the 210 ft. radio telescope at 
Parkes, Australia, by Hazard and associates (Reference 53) at 136, 410 and 1420 Mc/s. This led 
not only to the so far most accurate determination of the position, but also to a determination of 
the inner structure of the source. It consists of two components, designated A and B. B is as- 
sociated with the star-like object, A with the jet (Figure 14). The source B seems to consist of a 
bright but small core with a width of 0V5, producing about 80 percent of the radio emission, and a 
halo of circular or elliptical shape with a width of about 7”. The source A seems to have a similar 
structure, with a somewhat elliptical shape with the major axis in the direction towards com- 
ponent B. According to new eclipse observations, the difference in position between the star-like 
object and source B is less than 0”.5.* This is so extremely small that one cannot doubt the identi- 
fication, in particular since the source A coincides with the optical jet within the same accuracy. 

From the lunar occultation measurements. 

Hazard and collaborators were able to derive 
the relative intensity of the two components at 
410 and 1420 Mc/s and at 136 Mc/s. The ratio 
at the latter frequency, however, is highly un- 
certain. Using these data and the intensities 
measured by Mills at 85 Me (Reference 32) and 
by Dent and Haddock at 8000 Me (Reference 54) 
one can draw the radio spectrum of the entire 
source and the more uncertain spectra of the 
individual components (Figure 15). As Dent and 
Haddock pointed out, the spectrum of the com- 
ponent B resembles a thermal spectrum where- 
as A shows a typical spectrum of magneto- 
bremsstrahlung (synchrotron emission). 

If one, however, follows their idea to ex- 
plain the B spectrum by thermal emission, one 
arrives at highly unlikely results. Under the 
assumption of thermal emission, one can 
derive the temperature, electron density and 
total mass of the gaseous matter, if one knows 
the distance and the angular diameter of the 
source. The observed redshift, interpreted as 
purely cosmological, leads to a distance of 
about 5 x 10 8 pc. 

*A. Moffett, private communication. 



Figure 14— Optica! and radio features of 3C 273. The 
"star" designates the position of the optically bright 
object; the large hatched area, the faint jet. The circles 
(A, B) show the positions of the related radio sources 
each consisting of a luminous core (inner circle) and a 
halo of elliptical or spherical shape. The triangle gives 
the edges of the moon during the 3 occupations in 1962, 
determining the high accuracy of the radio positions 
(see Reference 53). 
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The diameter of the core of B which emits 
80 percent of the radiation was measured to 
be 0,5". This corresponds to a linear diameter 
of 1.3 kpc. From the shape of the spectrum, 
one easily sees that an optical thickness equal 
to one is reached at about 200 Me. Below this 
frequency, where the slope is 2, the radiation 
emerges from an optically thick layer. The 
flux of 6 x 10" 26 W/m 2 c/s at 100 Mc/s yields 
a temperature of 4x 10 9 °K. Furthermore 
r = 1 at 200 Me/ s yields an electron density 
of 10 5 cm" 3 and a total mass of 10 12 M 0 . 

If one assumes that the real diameter is 
only 1/10 of the observed upper limit, then 
the temperature would have to be increased 
to 4 X 10 11 , the electron density to 10 7 cm" 3 . 
The total mass would be 10 11 M 0 . 

All these numbers are so large that an 
interpretation of the radio spectrum of com- 
ponent B as thermal emission seems to be 
very unlikely. On the other hand, an inter- 
pretation such as magneto-bremsstrahlung requires either an envelope of ionized gas with suf- 
ficiently high electron densities roughly on the order of N e = 10 7 cm" 3 which absorbs the lower 
frequencies (below 100 Me) or that the energy spectrum of the relativistic electrons must have a 
cut-off towards the lower energies. For example, in a field of 10" 5 gauss, the cut-off would have 
to occur at about 2 x 10 9 ev. 

Immediately after Maarten Schmidt’s finding of the large red shift in 3C 273, Greenstein (Ref- 
erence 55) reexamined the spectrum of 3C 48. This spectrum does not show a Balmer series, but 
applying a redshift of 0.368 corresponding to an expansion velocity of 111,000 km/sec, he could 
identify the forbidden lines of O n, Ne HI and Ne V. The strongest line in the spectrum then ap- 
peared at 3830 A, which belongs to ionized Mg with an unshifted wavelength of 2800 A. 

Recently Maarten Schmidt found other star-like radio sources with still larger redshifts: 

3C 47 with a redshift of 0.425 corresponding to a nominal velocity of 127,000 km/sec and 3C 147 
with a redshift of 0.54 corresponding to 162,000 km/sec. According to Schmidt and Matthews 
(Reference 56) the spectrum of 3C 47 is similar to 3C 48, but the angular diameter of the radio 
source is much larger, about 1 minute of arc corresponding to 400 kpc. This size corresponds 
already to the largest double source radio galaxies. One might speculate whether the star-like 
radio sources are just another phase in the evolution of a radio galaxy. The redshift of 3C 147 is 
still somewhat tentative, since the two strong lines could also be explained with a redshift of 
0.229. But this is less likely. 



Figure 15— Radio spectrum of 3C 273. The black dots are 
measurements of both components. The squares give the 
spectrum of component A; the open circles of component 
B. The line through the open circles shows how one could 
fit a thermal spectrum to it, but this interpretation seems 
highly unrealistic (see text). 
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Figure 16 shows the schematic spectra of 
the 5 star-like radio sources with identified 
redshifted spectra (3C 273, 3C 48, 3C 47, 
3C 147 and 3C 286) together with the spectra 
of the famous radio galaxies Cyg A and 3C 295. 
The peculiarity of 3C 273 is that the forbidden 
O n line (3727 A) is so weak that it hardly can 
be seen on the spectra and therefore was not 
listed in Schmidt T s original paper. Very re- 
cently, Greenstein and Schmidt (Reference 57) 
published an excellent detailed analysis of the 
radio sources 3C 48 and 3C 273. 

A few words should be said about the pos- 
sibilities that the redshift is not purely cos- 
mological, but rather gravitational. In that 



case, the quasars could be peculiar stars 
within the Milky Way, rather than far away ex- 
tragalactic objects. But a close examination 
of the implications of gravitational redshift 
rules this possibility out. The gravitational 
redshift depends on the ratio of mass to 
radius. Greenstein derives for 3C 48 with 
Z = 0.367 a radius R (in km) = 4 M/M 0 and a 
radius R(in km) = 9 M/M 9 for 3C 273 with 
Z = 0.158 

A star with about 1 solar mass and a rea- 
sonable surface brightness and a diameter of 
10 km would have to be so close that its proper 


Figure 16— Schematic spectra of the quasars 3C 273 , 48, 
47, 147 and 286, together with Cyg A and 3C 295. The 
ordinate is the redshift A\/\. The scale is logarithmic 
in 1 + AAA (adapted from Reference 55). Forbidden 
lines are represented by black dots. The numbers under- 
neath the dots indicate the relative strength of the lines 
(number 1 being the strongest line in the spectrum). The 
redshift of 3C 147 is labeled as tentative, but the most 
likely between two possibilities (Reference 56). The 
redshift of 3C 286 is extremely uncertain, since there is 
only one line observed and no other evidence available. 
The identification with the Mg II 2800 A line was pro- 
posed by Shklovsky. On the right the following data are 
given: 

column 1: the observed redshift A\A , 
column 2: the nominal "recession speed" v/c obtained 
by inserting A X/k into the Doppler formula, 
column 3: the observed visual magnitude. 


motion should have been observed. Jefferys 
of Yale University, however, reported (Ref- 
erence 58) that there is no proper motion more than half the mean error of the determination. He 
used 14 plates from 1887 through 1963. The absolute proper motions he obtained are about 


0.001 ± 0.0025 sec of arc per year. This yields a lower limit for the distance of 2 x 10 4 pc; there- 
fore, the object is outside of the galaxy. 


Another piece of evidence in the same direction is provided by Williams of the University of 
California in Berkeley. At the Symposium on Gravitational Collapse in Dallas 1963 he reported 
(Reference 59) about observations of the 21 cm line absorption of the galactic hydrogen situated 
between the source 3C 273 and the Earth. Unfortunately the source is in the constellation Virgo 
at the large galactic latitude of 64°. Therefore the line of sight does not intersect dense hydrogen 
clouds. But it seems that all hydrogen, which can be observed in emission in the surrounding of 
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the direction toward the source, appears in absorption, when the telescope is pointed on the source. 
The effect, however, is rather marginal. 

A strong argument against gravitational redshift is the fact that there are strong forbidden 
lines which can only originate in highly dilute gases, but hardly in an atmosphere of a small object 
with extremely high density. A further argument is that the spectral lines are rather sharp (their 
widths correspond to about 50 km/s). If the lines, indeed, originate in an object with a 10 km 
diameter, the variation of the redshift within the layer, from which the lines are emitted, would 
cause much broader lines. 

Therefore one can use the redshift as cosmological distance indicator with good confidence. 
Using these distances, one finds that these star-like objects are by far the brightest objects in the 
universe. Table 5 gives the nominal speed, the observed visual magnitudes and the absolute magni- 
tudes based on a Hubble constant of 100 km sec" 1 Mpc -1 and on an assumed linear relation between 
redshift and distance. Furthermore two values for redshift-corrected luminosities are given. The 
first (- 2 A\A) is the necessary minimum correction, the second (- 5 A\A) applies for a spectral 
energy distribution in a normal galaxy. Since the star-like objects appear to have high intensities 
in the blue and UV range, the necessary correction would probably lie in between the two given 
values. For comparison the corresponding data for the two giant radio galaxies Cyg A and 3C 295 
are also given, and so are the absolute magnitudes of the brightest observed elliptical galaxy and 
our own galaxy. It is evident that the star-like objects are about 100 times brighter than normal 
giant galaxies. Typical values for the total emission in the optical and the radio range were already 
given in Table 4. 


Table 5 

Apparent and Absolute Magnitudes of 3 Quasi-stellar 
Radio Sources and Two Radio Galaxies.* 


Object 

cz 

10 3 km/s 

m v 

M 

V 

M 

V 

AX 

- 2 X 

AX. 

M v - 5y 

3C 273 

47.5 

12.6 

-25.6 


-25.9 

-26.4 

3C 48 

no 

16.2 

-24.5 


-25.2 

-26.3 

3C 47 

127 

-18 

-23 


-24 

-25 

Cyg A 

1G.8 

14.1 

-21.5 


-21.6 

-21.8 

3C 295 

138 

20.9 

-20.1 


-21.0 

-22.4 

Brightest elliptical galaxy: 


- 22.7 

Our galaxy 




-21.0 


*Two red shift corrections on the absolute magnitudes have been 
given, in between which the true values can be supposed. For 
comparison, the absolute magnitudes of the brightest elliptical 
galaxy and of our own galaxy are also given. 


In 1963 Matthews andSandage (Reference 60) 
discovered light variations in 3C 48 and Smith 
and Hoffleit in 3C 273 (Reference 61). In 
3C 273 there seem to exist small cyclic vari- 
ations with a period of about 13 years; flashes 
during a week or a month during which the 
object is up to 1 mag brighter and in 1929 a 
sharp drop of 0’.’4 occurred, after which the 
object restored to normal brightness in 1940. 
These variations impose a severe limit on the 
size in which they can occur. The 13-year 
cycle imposes no problem, if 3C 273 has only 
a diameter of about 7 light years according to 
the optical model by Greenstein and Schmidt 
(Reference 57), but the shorter variations 
might occur in only local regions of the object, 
in which case the energy content reaches very 
high values. Several processes have been 
suggested but even a supernova event would 
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not change the magnitude of the object appreciably, so that for the present this remains an open 
problem. 

The sudden drop in 1929 took place in 7 months and must have occurred through the whole 
object, thus limiting the size of the body, where the optical radiation originates, to .3 pc, and re- 
sulting in a very high electron density in that gas cloud (Reference 57). 

What can be said about the lifetimes of these radio sources? Unfortunately only crude times 
within several orders of magnitude can be obtained. The outer edge of the jet in 3C 273 is 1.5x 10 5 
light years away from the star-like object. Since it is likely that the jet originated from the star- 
like object, we obtain a lower limit for the lifetime of this object of 1.5 x 10 s years. 

Among the brightest radio sources in the 3C Catalogue, most of which (namely 25) are identi- 
fied, we find 4 star-like radio sources. This corresponds to 15 percent. One might conclude from 
this that they either are in fact more seldom or that their lifetimes are shorter by a factor of 5 
compared to the radio galaxies. This method, however, can give a hint only. The latest identifica- 
tions indicate a ratio of about 2:1 between radio galaxies and quasars. 

If we assume a lifetime of 10 7 years as a likely value, then the star-like objects have to pro- 
vide a total energy in the form of relativistic electrons of 3 x 10 5 8 erg. In the case of the largest 
radio galaxies, this value goes up even to 10 60 erg. Since the efficiency of producing relativistic 
electrons cannot be very large, the actual requirements for the energy supply are still larger by 
at least two orders of magnitude. If one considers that 2 x 10 62 erg corresponds to the rest-mass 
energy of 10 8 solar masses, one will be able to appreciate the difficulty in explaining the energy 
supply for these strange objects. 

Radio Sources and Clusters of Galaxies 

Recent investigations by Pilkington (1964) show that there exists a significant correlation in 
position between discrete radio sources and clusters. Using the new 178 Me Cambridge survey 
and Abells catalogue of rich clusters, taking an area containing 435 sources, he found that 22 
sources had positions within 0.4 times the radius of the cluster. Only 3 out of the 22 can be ex- 
pected to be due to chance. A similar investigation using the catalogue by Mills, Slee and Hill 
(Reference 32), yielded 21 more coincidences, within .6 times the cluster radius, of which 5 might 
be due to chance. The coincidences tend to be close to the center of the cluster, which is in 
agreement with the fact that many D galaxies, often being close to the centers of clusters, are 
strong radio sources. Another point of evidence that the coincidences are real is the fact that the 
apparent radio flux from sources coinciding with nearby galaxies is generally larger than from 
those coinciding with far away clusters. There was no evidence that the richness of the cluster 
affects the probability of strong radio emission. A quite different picture is presented when we 
use only the identified radio sources (Reference 34). The results then obtained are: 50 percent 
of the parent galaxies occur in clusters with richness ^0 and about 33 percent in clusters with 
richness >1, but only 20 percent appear not to be associated with clusters of galaxies; at most 
they might belong to very small groups. The quasars do not seem to belong to clusters, based on 
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the poor statistical evidence. The stronger sources tend to occur in clusters of rich- 
ness 2. 

This shows very clearly that one has to be very cautious in exercising statistics with these 
objects. On one hand many faint clusters are not included in Abell's list and therefore reduce 
Pilkington’s coincidences (Reference 62). On the other hand, using only identified sources certainly 
does not present a homogeneous sample. One may conclude however that beyond any doubt there 
is a correlation between radio sources and clusters, the degree of correlation has yet to be de- 
termined by avoiding selection effects as much as possible. 

The Radio Spectra 

We will discuss here the spectra of the radio galaxies and quasars. The radio spectra of 
nearly all radio sources are very similar in shape and suggest strongly a great conformity re- 
garding the processes governing the radiation of these sources in the radio frequencies. A source 
radiating by means of the synchrotron radiation process would have a spectrum represented by 
F v ~ v 0 = v K 0 ' , provided that the energy spectrum of the relativistic electrons has the form 

N e (E) dE = Ae 8 ° dE (see Equation 82). Hence plotting log F v versus log v we expect to find a 
straight line. This is essentially confirmed by the majority of the spectra up to now obtained 
(References 26b, 63 and 64). The average value of the spectral index is found to be 0.74. Depar- 
tures from the straight line can be explained in several ways, e.g.: 1) the electron energy spectrum 
is not a power spectrum, 2) the radio spectrum may be the superposition of several components. 

1) Because of energy losses due to synchrotron radiation, ionization and Coulomb- 
bremsstrahlung, the energy spectrum of the electrons may become quite different from the spectrum 
of the injected electrons (the source spectrum). As these losses vary differently with energy for 
the different processes, the spectrum will be determined by the process responsible for the 
largest losses. If synchrotron radiation dominates, the spectral index will tend to a-a 0 + 1/2 ; if 
ionization is the main process, a-*a 0 - 1/2, with all transitions possible (Reference 63). Hence in 
sources with very large absolute power and high temperature brightness, where we might expect 
relatively strong magnetic fields, the spectra will be considerably steeper than for less strong 
sources, where other processes may take over, which all tend to flatten the spectra. 

Five of the known quasars tend to have very steep spectra, though a is never larger than 1.25, 
agreeing with a = a 0 + 1/2 , in the extreme synchrotron radiation case. Furthermore, nearly all 
spectral indices lie within the a 0 + 0.5, a 0 - 0.5 range (i.e., between 0.24 and 1.24). This might 
be interpreted as a confirmation of the dependence of the spectral index on the dominating process 
for energy losses of the electrons (Reference 65). This outline holds only for a process of con- 
tinuous supply of electrons. If the supply of relativistic electrons ends, the energy spectrum of 
the electrons will show a band, which will move to lower energies with increasing time. This 
would result in a very steep part in the radio spectrum from a certain frequency onward (depend- 
ent upon the magnetic field and time); this frequency would then also decrease with time. How- 
ever no sudden steepening has been observed as yet. The currently accepted values for the inter- 
galactic magnetic fields and the electron lifetimes would give a cut-off frequency of 10 4 - 10 s Me, 
for which no data are available as yet. 
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Kellerman suggests that if the electrons would be supplied from the intergalactic space, it 
would explain the absence of a cut-off and the remarkable similarity of the energy spectra of the 
relativistic electron source for many electrons in so strongly different radio sources. 

2) One can easily see that for a complex source with unequal components with different spec- 
tral indices, the curvature of the integrated source will tend to be positive (compare portions A 
and B of Figure 14 of the radio spectrum of 3C 273) even if the curvatures of the components are 
negative (for low frequencies the steep component will dominate, for high frequencies the flat com- 
ponent). This is in good agreement with observations by Howard, Dennis, Maran and Aller (Ref- 
erence 64); they found: 

a) sources with positive curvature tend to be double, halo-core, or complex and of greater 
linear extent than sources with negative curvature. It is quite possible that all sources with posi- 
tive curvature have a complex structure, 

b) sources with negative curvatures tend to be smaller in linear size and in separation, and 
tend to have higher luminosities and surface brightnesses (Reference 66). Five quasars were 
found to have spectra with large negative curvature as hence could be expected. For the other 
quasars more data are needed. 

As only few sources with positive curvature were found, and the great majority of all sources is 
double or complex, the components are generally expected to have a similar physical structure, 
and to have the same origin. 

Size and Separation of Radio Sources and Their Components 

The radio features of about all identified strong sources show a double or more complex 
structure. Maltby and Moffet (Reference 46) studied 174 extragalactic sources; they were able to 
resolve 75: 

13 did not have a complex structure, 

15 had two equal components, 

40 had two or more components of unequal intensity, 

7 had a bright core source and a faint halo. 

Correcting for the spatial distribution, which makes us see several double- sources edge-on, we 
find that 80 - 85 percent of all resolved sources have double features, 10 percent are of the halo- 
core type, and only 5 to 10 percent are single. These percentages hold also for identified sources 
(Reference 34). 

The double structure of nearly all strong sources is centered about the nucleus of the parent 
galaxy, but the components do have separations up to 360 kpc (Centaurus A). Maltby, Matthews and 
Moffet (Reference 67) suggested that due to some catastrophe occurring in the nucleus of the galaxy, 
a radio source of small diameter (1 kpc) had been ejected from the galaxy nucleus. While moving 
away from the galaxy the source might then expand to more than 100 kpc when having a distance of 
150 kpc to the nucleus. 
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Many of these galaxies have a rotational symmetry, recognizable by ring- like absorbing dust 
bands, and sometimes actual rotation has been observed (Reference 39). Where detailed optical 
observations were possible faint extensions in the direction of the axis of rotation have been found. 
These are very peculiar features in a rotating body, and can only be explained by the presence of 
a magnetic field aligned with the axis of rotation, and along which a plasma has been ejected. In 
many cases the main axis coincides approximately with the rotational axis of the galaxy and the 
components are extended in the direction of that axis. 

Sometimes the components are connected by an emissive bridge and show rather sharp ex- 
ternal edges. This might be associated with the shock front of a plasma moving through a 
magnetic field, away from the galactic plane (Reference 68). 

Since it seems obvious that the magnetic field structure of the galaxies is decisive for the 
structure of the sources, many experiments have recently been carried out to detect any signifi- 
cant polarization in the radio region. As the main emissive process in these sources is believed 
to be synchrotron radiation, the polarization should give some information regarding the strength 
and the direction of magnetic fields. 

Optical polarization experiments have been carried out by Baade on Messier 87 (Refer- 
ence 49) radio experiments by Westerhout (Reference 69) Bracewell and associates (Reference 
40) Mayer and associates (Reference 48), Morris and associates (References 70, 71 and 72). 

In nearly all studied cases linear polarization of both components was found. The difference 
between the angles of polarization of both components was generally smaller for sources with 
small separations, than for sources with larger separations. This is in agreement with the 
model of ejected matter from the galaxy: we expect the two components to behave rather 
independently from each other when they have moved a considerable distance from the galaxy; 
close to the galaxy there will probably still be a large interaction between the two components 
(References 70, 71 and 72). There seems to be no simple relation between the angle of 
polarization and the position angle of the main axis of the double source, although there is a 
tendency for a 90° difference. This would then give the result predicted by the model: an 

alignment of the magnetic field along the axis of the source, and approximately along the axis 
of rotation of the galaxy. 

3C 273 has two radio components, one halo-core type connected with the stellar-like object, 
the other component relates to the optically detected jet (see Figure 14). The source 3C 48 is 
single. It has a very small angular diameter (<1 ,T ). 3C 286 and 3C 196 are both probably halo-core 
sources, 3C 147 may be a halo-core or a double source (Reference 56). Hence, just like in the 
case of the optical spectra, the radio structures do not show any uniformity among the quasars. 

The linear sizes of the quasars and their components are generally much smaller than the other 
strong extragalactic radio sources, as the N galaxies and the D galaxies, which size up to 80 kpc 
in their optical dimensions. The diameter of the radio components ranges from <3 kpc for 3C 48 
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to 10 kpc for the halo of the 3C 273 B component. The separation of the components of 3 C 273 is 
about 40 kpc. 3C 47, the quasar with the smallest luminosity in both the optical and the radio range 
up to now detected, has a distance between the components of 250 kpc. This is a distance which is 
already in the same order of magnitude as the separations between radio components of large radio 
galaxies. 

Cosmological Implications of Radio Source Counts 

The accurate determinations of flux densities at a frequency of 178 Me by Scott, Ryle and 
Hewish (Reference 73) have confirmed the earlier findings that the relationship between the number 
of radio sources and their flux densities is steeper than one would expect from cosmological 
models, which are based on the assumption of homogeneous distribution of radio galaxies in space 
without evolutionary effects. 

Scott and Ryle (Reference 74) have shown that the number of sources plotted versus their 
observed flux in a double logarithmic scale can be represented by a straight line with a slope of 
1.8.* This holds true for a flux range from 2 to 100 x 10" 26 Watt m‘ 2 (c/s)" 1 at 178 Me. For 
smaller fluxes below 10“ 26 W m" 2 (c/s)" 1 a flattening of the slope is indicated in an investigation 
by Hewish (Reference 75). He applied a special method proposed by Scheuer which allows the 
extension of the search below the limit of the observation of individual sources. The flattening of 
the slope was confirmed in direct observations with an instrument of greater resolution for a 
small area of the sky (Reference 76). 

A slope of 1.8 in the number-flux relation cannot be explained with the models derived from 
relativistic cosmology or steady-state cosmology (References 77, 78 and 79) without introducing 
evolutionary effects. There is only little quantitative difference between the number-flux re- 
lations calculated for various simple cosmological models, but they all show a curved relation 
which starts at large flux densities with a slope of 1.5 and then bends over to smaller values. 

In these investigations the distribution of radio luminosity can reasonably be represented by 
an effective luminosity. This takes into account that the brighter sources are sampled from a 
larger volume if they are counted as function of the flux density. The luminosity functions derived 
by Ryle show that the effective total emission in the radio frequencies (10 - 10 4 Me) is in the range of 
10 43 to 10 44 erg sec" 1 . If there is a high percentage of quasars among the observed radio sources 
this value for the effective luminosity might somewhat increase. 

The requirements for evolutionary effects in the radio sources have been investigated by Oort 
(Reference 80) (1961) and Davidson (Reference 79). In order to represent the observed slope of 1.8 
one can either assume that radio sources were in the past more numerous per unit volume of the 
universe or that their effective luminosity was considerably larger in the early times of the evolu- 
tion of the universe. Of course both possibilities can contribute in part. 


•In order to have a positive slope it is assumed that the flux scale is increasing to the left side in the number-flux diagram. 
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From the formulas given by Priester (Reference 78) one can easily derive a lower limit for 
the requirements of the evolutionary effect if one requests obtaining a slope steeper than 1.5. We 
shall give here the relation for the Einstein-de Sitter model with euclidean metric and a cosmo- 
logical constant A = 0. In this model the expansion is proportional to (t 0 -t) 2/3 where t 0 is the 
present age of the universe as defined in this model and T the time counted from the present time 
into the past. T is used to characterize the travel time for the radio waves from the source until 
they reach the observer. 

The requirement for the time dependence of the number densities of radio galaxies per unit 
volume is 



where the exponent a > 4. For the luminosities the requirement is less severe 



( 102 ) 


(103) 


where b > 8/3, n 0 and P 0 are the present values for the number density and luminosity derived from 
the closest objects. 

If we assume an age of the universe of 10 10 years we find that for times of 3, 5 or 7 x 10 9 
years ago, the required number density of radio galaxies was more than 4, 16 or 120 times larger 
than at present. The corresponding values for the required luminosities are 2.6, 6.3 or 30 times 
the present luminosity for 3, 5 or 7 x 10 9 years ago, respectively. Because of the singularity of 
the model for T ~ t 0 one is, of course, not permitted to extend this result up to that time. From 
the numbers given it is apparent that in particular for the luminosities, the required increase in 
the past does not seem to be unreasonable, in particular since there is also the possibility that 
both the numbers and the luminosities were larger then. This, of course, would lessen the re- 
quirements on each of them correspondingly. In any case it must be kept in mind that the estimates 
given above can only provide some insight into the order of magnitude required for evolutionary 
effects. 

Physical Processes 

We will discuss here only the strong sources, following the discussion on this subject as pre- 
sented by M. Ryle at the URSI Assembly 1963. We have to divide the processes into two classes, 
relating to different problems: 1) source of energy, necessary to produce these high luminosities, 
2) processes to explain the spatial distribution of the radio emission or of the magnetic fields. 
Further the emitted power and the shape of the spectrum must be explained if one assumes the 
existence of a source of energetic particles and a magnetic field. 
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1) The main problem with the strong sources is that they require, even with modest estimates 
of their lifetimes, very large total energies. Theories which try to provide the energy from ex- 
ternal sources, for instance from collisions between galaxies or even with large clouds of anti- 
matter, are less favoured nowadays than those based on processes within the galaxy. All have to 
explain a total energy of the order of 10 58 - 10 60 ergs, supplied in the form of relativistic electrons. 
This value is obtained as follows: the lifetime of the relativistic electrons causing the radio fre- 
quency radiation is of the order of 10 7 years, if the magnetic field has a field strength of H = 10~ 5 
gauss and if the initial energy of the relativistic electrons is of the order of 10 9 ev. As these 
sources emit between 10 41 and 10 45 erg sec" 1 in the radio range and up to 10 46 erg sec" 1 in the 
optical range, the energy supplied in the form of relativistic electrons is of the order of 10 5 9 erg, 
if the lifetime is correctly estimated. From this one can conclude that the total energy for the 
largest sources is likely to exceed 10 61 erg since there will surely be only a small percentage of 
the total energy stored in relativistic electrons. 

Supernova explosions cannot explain this amount of high energy electrons, if they occur at the 
same rate as they explode in our galaxy, which is 1 every 300 years. One explosion yields about 
10 51 erg, of which only a very small portion could be converted into high energy electrons. In 
10 7 years this would only yield 10 56 erg. 

A collision between two giant galaxies would not be able to supply the required energy either. 
Let their masses be 10 11 solar masses, and their relative velocity 500 km/sec. Then the collision 
energy is 6 x 10 59 erg but here too only a very small part of the energy would be transfered into 
high energy electrons. 

Proposed mechanisms for radio galaxies and quasars involve supernovae explosions or the 
gravitational collapse of a superstar. 

Shklovksy (Reference 9a and 9b), Burbidge (Reference 81) and Cameron (Reference 82) argue 
convincingly that a plausible way to obtain a sufficient number of relativistic electrons is to assume 
an excessively higher supernova rate during some periods within the evolution of the galaxy. 
Shklovksy assigns this to an early stage of the galactic evolution. Burbidge assumes that a chain 
reaction of Type I supernovae in the very dense central regions of a galaxy can provide the neces- 
sary frequency of explosions. In this close-packed cluster of old stars the explosions are trig- 
gered by the shock waves from the exploding neighbor stars. It seems not implausible that the 
required very high star densities might be found in the nuclei of D galaxies, N galaxies or perhaps 
in quasars. 

Strong radio galaxies are usually of the D and N types, which can be considered as a special 
form of elliptical galaxies. This fact favors also the suggestion by Cameron (Reference 82) who 
assumes that intense supernova activity can arise in these galaxies during the star formation 
process in a large gas cloud. The more massive stars contract quickly and all spend about 3 x 10 6 
years on the main sequence in the Hertzsprung-Russell diagram. Subsequent evolution will be 
rapid and will lead to Type II supernovae within a relatively short period (of the order of 10 7 
years). 
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Hoyle and Fowler (Reference 83) try to derive the required energy from a gravitational col- 
lapse of an excessively massive superstar (10 4 5 6 7 to 10 8 * 10 11 solar masses). This suggestion is in par- 
ticular of importance in connection with the quasars. The difficulties of this idea are that the 
contraction will lead to a rotational crisis, a faster and faster rotation which will break up the 
object into smaller objects. On the other hand this provides a possibility for explaining the varia- 
tions in the optical emissions, as Hoyle argued at the Symposium on Gravitational Collapse in 
Dallas 1963. If there are large numbers of subcondensations formed, which will reach their 
luminous stage statistically distributed, one might expect noise-like light variations. A more 
detailed account on the ideas of explaining the physics of radio galaxies or quasars might be found 
in the reports on the Dallas symposium by Chiu (Reference 1) and Weidemann (Reference 84) or 
in the proceedings of the symposium edited by Robinson, Schild and Schucking (Reference 85). 
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